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EXCLUSION PROCESSES IN HIGHER DIMENSIONS: 
STATIONARY MEASURES AND CONVERGENCE 

By M. Bramson^ and T. M. Liggett^ 

University of Minnesota and University of California, Los Angeles 

There has been significant progress recently in our understand- 
ing of the stationary measures of the exclusion process on Z. The 
corresponding situation in higher dimensions remains largely a mys- 
tery. In this paper we give necessary and sufficient conditions for a 
product measure to be stationary for the exclusion process on an 
arbitrary set, and apply this result to find examples on Z'^ and on 
homogeneous trees in which product measures are stationary even 
when they are neither homogeneous nor reversible. We then begin 
the task of narrowing down the possibilities for existence of other 
stationary measures for the process on Z''. In particular, we study 
stationary measures that are invariant under translations in all direc- 
tions orthogonal to a fixed nonzero vector. We then prove a number 
of convergence results as t ^ oo for the measure of the exclusion 
process. Under appropriate initial conditions, we show convergence 
of such measures to the above stationary measures. We also employ 
hydrodynamics to provide further examples of convergence. 

1. Introduction. In this paper we consider the stationary measures and 
convergence for the exclusion process rjt on the countable set S with (stochas- 
tic) transition probabilities p{x,y). This is the Markov process on {0, 1}"^ 
in which a particle at x £ S attempts to move to y £ S at rate p{x, y). The 
move takes place if y is vacant and does not take place if y is occupied. 
The exclusion process has been very heavily studied since its introduction 
by Spitzer [18] in 1970. For an account of many known results about it, see 
the books by Liggett [13, 14], and Kipnis and Landim [8], and the references 
therein. 

In spite of the tremendous success that the study of the exclusion pro- 
cess has seen, relatively little was known about the existence and structure 



Received March 2004; revised October 2004. 
^Supported in part by NSF Grant DMS-02-26245. 
^Supported in part by NSF Grant DMS-03-01795. 
AMS 2000 subject classification. 60K35. 

Key words and phrases. Exclusion process, stationary measures, hydrodynamics. 

This is an electronic reprint of the original article published by the 
Institute of Mathematical Statistics in The Annals of Probability, 
2005, Vol. 33, No. 6, 2255-2313. This reprint differs from the original in 
pagination and typographic detail. 



1 



2 



M. BRAMSON AND T. M. LIGGETT 



of nonexchangeable and nonreversible stationary measures for the transla- 
tion invariant exclusion process on Z'^ until 2001. (The only truly general 
result that was known is that all extremal translation invariant stationary 
measures are of product form — see Theorem 3.9 of Chapter VIII in [13].) 
This is still the case if d > 1. In one dimension, however, the papers by 
Ferrari, Lebowitz and Speer [5], Bramson and Mountford [2] and Bramson, 
Liggett and Mountford [1] made significant progress on this problem. The 
present paper is a first attempt to develop a corresponding theory in higher 
dimensions. 

We begin by summarizing some of the one-dimensional results, since they 
provide a model for what we would like to prove in higher dimensions. The 
main assumptions are S = Z, p{x,y) = p{y — x), J2x \^\pi^) < °*^) the 
process is irreducible. [Throughout this paper, irreducibility will mean that, 
for each x,y € S, there is positive probability that the Markov chain with 
transition probabilities p{-,-) goes from x to y in some number of steps.] 
First, two definitions: A probability measure fi on {0, 1}^ is said to be a 
blocking measure if it concentrates on configurations rj satisfying 

^ r/(a;) < oo and [1 — ri{x)] < oo, 
x<0 x>0 

and it is said to be a profile measure if 

lim /i{?7 : ry(x) = 1} = and lim fi{rj -.rj^x) = 1} = 1. 

x— >— oo x— »+oo 

Every blocking measure is a profile measure, but not conversely. As usual, 
if < /3 < 1 , will denote the homogeneous product measure with density 
P- 

Without loss of generality, we may assume J2x xp{x) > 0. We then have 
the following results: 

1. Either the extremal stationary measures consist exactly of 

(i) K,0<p<l} 
or of 

(ii) {i/p,0<p<l}U{/x„,nGZ}, 

where /io is a profile measure, and /i^ is the shift of /io by n units. 

2. If X]a;3;p(x) =0, then case (i) occurs. 

3. If ^x^Pi.^) > ^i^d p(-) has finite range, then case (ii) occurs and /io is 
a blocking measure. 

4. If > 0? P(^) ^'^d 'P{~x) are decreasing for x = 1,2, . . . ; p[x) > 
p{—x) for all X > 1; and Z]x<o^^P(^) < then case (ii) occurs and hq 
is a blocking measure. 

5. If Yl,x<ox'^P{x) — there exists no stationary blocking measure. 
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Statements 1, 4 and 5 were proved in [1], statement 2 appears on page 391 
in [13], and statement 3 is due to [2]. An important open problem that 
remains in one dimension is to determine whether nonblocking stationary 
profile measures ever exist. 

While the above description is the starting point for our study of station- 
ary measures in higher dimensions, one difference between one and higher 
dimensions can be seen immediately in case (ii) above. In one dimension, the 
extremal spatially inhomogeneous stationary measures are naturally indexed 
by a discrete parameter. In higher dimensions, the natural parametrization 
is continuous, rather than discrete. Here is a simple example in two di- 
mensions that illustrates this point. Suppose p{ei) =p{e2) =p(— 62) = 1/3, 
where ei = (1, 0) and 62 = (0, 1). Then, for each s G i?, an extremal station- 
ary measure iXg is obtained in the following way. Letting x = {x^^\x^'^'^) and 
s = n + a for an integer n and < o" < 1, //g is the measure for which r/(x) = 1 
if x^^^ > n, ri{x) = if x^^^ < n, and 'r]{x) is i.i.d. with density a if x^^^ = n. 

We turn now to a description of the results of the current paper. Pre- 
cise statements and definitions appear later in the relevant sections. For 
the moment, S is general. The first problem is to determine exactly when 
product measures are stationary for the process. Suppose that a{x) satisfies 
< a{x) < 1 for all x £ S, set 

7r(x) = a(x)/(l — a{x)), 

and let Uq. be the product measure on {0, l}'^ with marginals Ua{r] : 'i]{x) = 
1} = a{x). It is not hard to check that z/q is reversible for 7]t if and only if 
7r{x)p{x,y) = Ti{y)p{y,x) for all x,y G S — see page 34 in [6], for example. 

Here we are interested in the issue of stationarity of z^q,, rather than re- 
versibility. Theorem 2.1 on page 380 in [13] says that Va is stationary for 
rjt if either (a) p{x,y) is doubly stochastic and a{x) is constant, or (b) 
7r{x)p{x,y) = Tr{y)p{y,x) for all x,y £ S. As mentioned above, in the latter 
case, Va is in fact reversible for the process. In the former case, is re- 
versible if and only if p(x,y) is symmetric. In Section 2 we will show that an 
appropriate combination of conditions (a) and (b) is necessary and sufficient 
for Ua to be stationary for rjt. 

Section 3 applies this necessary and sufficient condition to the case of a 
translation invariant system on S = Z'^. Under a mild assumption. Theorems 
2 and 3 combine to show that z/q is stationary if and only if there is a 
V £ so that 7t{x) = 7r(0)e^^'^^ and p{z) = e^'^''"'^p{—z) for all z such that 
{z,v) 7^ 0. This allows us to construct many examples of stationary product 
measures that are neither homogeneous nor reversible. (Homogeneity of Va 
corresponds to a being constant.) 

In Section 4 we begin by considering several examples in which p{x,y) 
are the transition probabilities for a random walk on a homogeneous tree S. 



4 



M. BRAMSON AND T. M. LIGGETT 



Again, using the results of Section 2, we find examples of stationary product 
measures that are neither homogeneous nor reversible. Then, we describe an 
example on a rooted tree that exhibits a type of phase transition, where the 
extremality of the stationary product measures among stationary measures 
depends on the value of a parameter. 

In Sections 5 and 6 we return to S = Z'^ and study stationary measures 
that are invariant under translations in all directions orthogonal to a fixed 
nonzero v Z'^. We suspect that all extremal stationary measures have this 
property for some v G R'^, but are quite far from being able to prove this. A 
measure that has this partial homogeneity property and has density tending 
to in the —v direction and to 1 in the +v direction will be called a v-profile 
measure. Theorem 4 in Section 5 asserts that there are no stationary f-profile 
measures unless v and the drift vector xp{x) have nonnegative inner 
product. (Under an irreducibility assumption, "nonnegative" can be replaced 
by "positive.") Other approaches to this and related results under varying 
assumptions appear in Theorem 5 of Section 5, Corollary 2 of Section 6 and 
Theorem 10 of Section 8 (the last based on hydrodynamics). This rules out 
half of the possible directions v in which stationary u-profile measures can 
exist. 

Theorem 6 in Section 6 provides a way of showing that, if there exists a 
continuous family of extremal stationary u-profile measures, then there can 
be no other extremal stationary t!-profile measures. This is applied in Corol- 
lary 1 to examples from Section 3 in which families of stationary product 
measures are known to exist. Theorem 7 of that section provides some rea- 
sonable assumptions under which any extremal stationary u-profile measure 
is a blocking measure, in the sense that it concentrates on configurations tj 
such that 



Theorem 8 of Section 7 provides sufficient conditions for convergence to a 
mixture of stationary f-profile measures when the initial measure is u-profile. 
The main assumption (besides the existence of an appropriate one-parameter 
family of stationary u-profile measures) is that the initial measure have a 
finite expected number of I's in the "negative half" of each strip oriented in 
the V direction (with bounded cross section in the perpendicular directions), 
and a finite expected number of O's in the positive half of that strip. 

Section 8 presents some applications of hydrodynamical results of Re- 
zakhanlou [16] to obtain explicit convergence results when the initial set of 
I's is a wedge in Z^. In this case, the initial distribution is not ?;-profile. 
Another example with such an initial distribution is discussed in Section 9. 
It provides an illustration of how the asymptotic behavior of the process can 
depend on more than just the mean of its transition probabilities. 
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This paper provides only an initial understanding of the stationary mea- 
sures and of convergence behavior for the exclusion process in more than 
one dimension. There are a number of problems in this context that we are 
not able to solve at this time. These are discussed in Section 10. There, 
we will also compare our results and speculations with the much easier and 
well- understood situation for independent particle systems on Z'^. 



2. Product stationary measures: the general case. In this section both 
S and p{x,y) are general. Our objective is to give a necessary and sufficient 
condition for a product measure to be stationary for the exclusion process. 
This is done in the following theorem. 

Theorem 1. Let ai,a2, ■ ■ ■ be the distinct values ofa{x), 7rj = aj/(l — 
aj), and Ci, C2, . . . he the partition of S determined by Ci = {x £ S : a{x) = 
Oj}. The following is a necessary and sufficient condition for to he sta- 
tionary for rjt : 

(a) T.y&c, P{x^ y) = T>y€C, P{y^ ^) for all xeCi and all i, 
and 

(b) 7rip{x, y) = TTjp{y, x) for all x £ Ci,y G Cj and all i^ j. 

Proof. Following the proof of Theorem 2.1 on page 380 in [13] up to 
display (2.7), we see that z/q is stationary if and only if the right-hand side 
of (2.7) is zero for every finite Ac S. Dividing by the product YlueA^i'^)^ 
we see that a necessary and sufficient condition for stationarity is that 

V- a(j/)[l - a{x)]p{y, x) - a{x)[l - aiy)]p{x, v) ^ 
^ a(x) 

for every finite AdS. Now fix a u G S" and consider (1) with A = {u}. This 
becomes 

(2) [1 - a{u)] ^ a{y)p{y, u) - a(u) ^[1 - a{y)]p{u, y) = 0. 

y y 

(Note that the terms corresponding to y = n in the above sums, which in 
principle should not be included, cancel out.) Next, fix distinct v,w £ S and 
consider (1) with A = {v,w}. This becomes 

a{v) 

(3) 

^ q(?/)[1 - (^W)\p{y,'w) - a{w)\y - a{y)]p{w,y) _ ^ 
^ a(w) 

y^v,w ^ ' 
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Subtracting (2) with u = v and (2) with u = w from (3) [after dividing these 
expressions by a{v) and a{w) resp.], (3) becomes 

(4) [a{w) — a{v)] [■k{w)p{w, v) — 'k{v)p{v, w)] = 0. 
Therefore, (1) holds for all singleton and doubleton A^s if and only if 

(5) (2) holds for all u and (4) holds for all v,w. 

Next, we will check that (5) implies (1) for all A, which gives us the 
equivalence of (1) and (5). Let A be any finite subset of S. By (2), the 
left-hand side of (1) is 

^ a{y) [1 - a{x)]p{y, x) - a{x) [1 - a{y)]p{x, y) 

a(x) 

which can be written as 

. N7r(y)p(?/,2;) -7r(x)p(a;,y) 

Interchanging the roles of x and y above and adding the two expressions, 
we see that the left-hand side of (1) is 

which is by (4). 

Finally, we note that (5) is equivalent to (a) and (b) in the statement of 
the theorem. Statement (b) in the theorem is equivalent to (4). Given this, 
it suffices to show that statement (a) is equivalent to (2) for all u. But the 
left-hand side of (2) for uG d can be written as 

{1 - ai)^{l - Oj) [TTjp{y,u) -TTip{u,y)]. 

The summands above for j ^ i vanish by (b), so (2) for this u is equivalent 
to 

i'^iPiy^ - T^iPiu, y)] = 0. 

But this is statement (a), so the theorem is proved. □ 

3. Product stationary measures on S = Z'^. In this section we assume 
that S = Z'^ and p{x,y) = p{y — x) for all x,y. We will show later that, 
under a minimal assumption, the vr corresponding to a stationary product 
measure must have an exponential form. [Recall 7r{x) = a{x)/{l — a{x)) if 
the product measure is z^q.] In the next result, we obtain an easily checked 
necessary and sufficient condition for stationarity in this case. 
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Theorem 2. Suppose that Tr{x) = 7r(0)e^^''"^ for all x and some v G R'^. 
Then, is stationary for rjt if and only if 

(6) p{z) = e^'^'^'^pi-z) for all z such that {z, v) / 0. 

Proof. Condition (a) of Theorem 1 says in this case that 

^ [p{x,y) -p{y,x)]={). 

y- {y,'")={x,v) 

Letting z = y — x and using translation invariance, this becomes 

^ \p[z)-p{-z)\=^, 

z : {z,i;)=0 

which is automatically true. Condition (b) of Theorem 1 says that 

e(^:''> p(y -x) = e^y^'"^ p{x - y) whenever e<^'^> / e^^^'"^ . 

Again, letting z = y — x,we see that, in this case, (6) is equivalent to (b) of 
Theorem 1. □ 

Remark. Theorem 2 makes it easy to construct examples of exclusion 
processes with a significant number of nontranslation invariant stationary 
measures of product form, including many that are not reversible. For ex- 
ample, suppose d = 2, 

P(.ei)=pi, p{e2)=P2, p{-ei) = qi, p{-e2) = q2 

and p{z) = for all other z's. Here ei,e2 are the standard basis elements in 
Z'^ . Assume that pi ^ qi and P2 7^ 92 • Then the stationary product measures 
with density 1/2 at the origin are exactly the four corresponding to the 
following: 

7r(x) = l, 7r(2;) = 

<x) = ip2/q2f"' , vr(x) = ipi/qif' {p2/q2f"' , 
where we have written x = {x^^\x^'^^). Only the last of these is reversible. 

Next, we will show that, under a weak assumption (that is trivially nec- 
essary), vr must have the exponential form assumed in Theorem 2. Most of 
the work is contained in the proof of the following statement. 

Proposition 1. Suppose a is such that is stationary for rjt, and 
u ^ Z'^ satisfies p{u) > 0. 

(a) If p{—u) = 0, then tt{x + u) = tt{x) for all x £ Z'^. 



8 



M. BRAMSON AND T. M. LIGGETT 



(b) If p{—u) > 0, then for each x G Z"^, either 

7t{x + nu) = 7r{x) for all n£ Z, 



or 



7r{x + nu) = 7r{x)X 



n 



for all n(z Z, 



where 



(7) 



A 



p{-u)' 



Proof. Suppose first that p{—u) = 0, and fix an x. Then the equahty 
in (b) of Theorem 1 cannot hold with y = x + u since one side is zero and 
the other is not. Therefore, ■k[x) = 7r(x + u). This proves part (a) of the 
proposition. 

The proof of part (b) is significantly harder. Suppose p{—u) > 0, and 
define A as in (7). By Theorem 1(b), 

(8) for each x G Z*^, 7r(x + u)/7r(x) = 1 or A. 

Therefore, for each x and n, tt{x + nu)/TT{x) is a power of A, and the sign 
of the power is the same as the sign of n. Our problem is to show that, for 
each X, the power is identically n or identically 0. To do so, we will need 
to use part (a) of Theorem 1 as well. Note that the result is immediate if 
A = 1, so we may assume that A 7^ 1. By interchanging the roles of u and —u 
if necessary, we may assume that A > 1. 

As mentioned above, we will need to use the property in Theorem 1(a). 
This involves sums in which p{x,y) and p{y,x) both appear for x,y with 
the same value of vr. It will be important to relate these two transition 
probabilities so that one sum can be related to the other. That is the purpose 
of the next observations. We will call the two sites of interest z and w instead 
of X and y to emphasize that, for the moment, they are fixed, rather than 
being summed over. In (10) and (11) below, we note that, depending on how 
z and w are placed relative to the values of vr at these and "neighboring" 
points, one of the transition probabilities, p{z,w) and p(w,z), is A times the 
other. 

Consider two distinct points z,w Z'^ for which '^{z) = 7r{w). If p{z, w) > 
and p{w, z) = or vice versa, then by part (a) of this proposition, 

(9) Tr{z + nu) = 'k{w + nu) for all n£ Z. 

Suppose now that p{z,w) > and p{w,z) > 0. If Tr(w + u) ^ T^iw), then 



by (8) 



tt{w + u)/tt{w) = A. 
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If also 'ir{z + u) = 'k{z), then ^{w + u) ^ -k{z + u), so by Theorem 1(b), 

p{z, w) = Xp{w, z). 
Similarly, if tt^w — n) / tt^w) and it{z — u) = vr(z), then 

p{w, z) = Xp{z, w). 
Now let c be a value of the function 7r(-), and write 

A = {z: Tr{z — u) = A~"^c, n{z) = c, 7r(z + u) = Ac}, 
B = {z: tt{z — u) = X~^c, 7r{z) = c, 7r{z + u) = c} 

and 

C = {z: it{z — u) = c, 7r{z) = c}. 

Note that by (8), 

AuBuC = {z:tt{z) = c}. 

li wGAUB and z £C, (9) is false (for n = —1), so p{z,w) and p{w,z) are 
either both positive or both zero. If they are both positive, then we know 
from the previous paragraph that p{w,z) = Xp{z,w) and, of course, this is 
also true if they are both zero. So, it is always the case that 

(10) w£AuB,z£C implies p{w, z) = Xp{z,w). 
Similarly, 

(11) w £ B,z £ A implies p{w, z) = Xp{z,w). 
By Theorem 1(a), since {z : tt{z) = c} = ALIBLIC, 

(12) \pi^^'^) -p(.w^z)] + Y\pi^^w) -Pi'^^^)]=^^ weB. 

zeAuc z&B 
Combining (10), (11) and (12) gives 

(13) (A-1) P(.w,z) = X^iPi^^w) -p{w,z)], weB. 

z£AUC z€B 

Note that w £ B implies z = w + u £ C, so that (13) implies 

(14) ^[p(z,ii;)-p(u;,z)]>(l-A-i)p(n)>0, w£B. 

z&B 

If d = 1, we can now complete the proof easily, since in this case, B is 
finite by (8). The sum on u; G i? of the left-hand side of (14) is zero, so the 
left-hand side itself is zero as well. It follows that B = 0. Similarly, 

{z : 7r{z — u) = c, 7r{z) = c, Tr(z -|- u) = Ac} = 0. 
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Together with (8), this completes the proof of the second part of the propo- 
sition in this case. 

An extra argument is needed in higher dimensions. The objective is still 
to show that i? = as a consequence of (14). Once this is done, the proof 
of the proposition is completed as in the previous paragraph. 

To do so, assume B ^ and let D C B be finite and nonempty. Sum (14) 
for w £ D. The part of the resulting sum on the left that corresponds to 
w,z £ D vanishes. Therefore, making the change of variable z = w + x and 
dividing by the cardinality of D, one obtains 

#{weD:w + xeB\Dy 



(l-A-i)p(n)<^ 



\p{-x) - p{x)\- 



\D\ 

Now let F C Z'^ he finite. It follows from the above inequality that 

n ( \^\^\ f ^ ( M . i^{{x.w)£F ^D:w^x£B\D} 
(1-A )VKn)<\^\v{-x)-p{x)\^ — . 

To get a contradiction, it will suffice to find, for any finite F, an increasing 
sequence of finite, nonempty subsets of B so that 

(15) ,^,^,*{M^F.B n,^x£B\D^^ 

" l^nl 

Once this is done, it will follow that the sums on the right-hand sides of 
the previous two displays can be made arbitrarily small by choosing F suf- 
ficiently large. 

Let -Do = 0) -Ci be any finite nonempty subset of and define recursively 
Dn+i = {w ^ x:x £F^w £ Dn, w + x £ B} U Dn- 
If (15) fails for this sequence, then there is an e > so that 

(16) #{{x,w)£FxDn:w + x£B\Dn}>e\Dn\, n>l. 
To show that this is not possible, consider the mapping 

{x, w) ^ X + w. 

This maps {{x, w) € F x Dn :w + x £ B\ Dn} into \ Dn, and is at most 

|F| to 1. Therefore, by (16), 

(17) e|D„|<|F||Z)„+i\Z),|, 

so that \Dn\ grows exponentially rapidly. On the other hand, Dn is contained 
in the union of |-Di| balls of radius Kn, where K = max{||j;||,x G F}, so 
\Dn\ can grow at most like K'n'^ for some constant K'. This contradicts the 
exponential growth, so (15) holds. □ 

Next we extend the statement of Proposition 1 to several values of u. 
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Proposition 2. Suppose a is such that is stationary for rjt , and that 
p{ui) > for 1 <i < k. If p[—Ui) > 0, set 

(18) A. - ^("^^ 



P{-Ui) 

while if p{—Ui) = 0, set Aj = 1. Then 

(19) 7r(niniH h n^-Ufc) = 7r(0)/)"^ • • • 

where for each i, pi = \i or 1. 

Proof. We prove this for k = 2, since the proof for general k follows 
by iterating the argument. So, suppose p{ui) > and p{u2) > 0. By Propo- 
sition 1, for each x, 

' TT{x + mui){pmT, 



TT(x + mui+nu2) — ^ , , N/ 

[7r{x + nu2)[a.n 

where pm = A2 or 1 for each m, and cj„, = Ai or 1 for each n. Putting n = 
and then m = in the above expression, and using the resulting expressions, 
we see that 

7r(x + mui + nu2) = 7T{x){ao)'^{pmT = T^{x){po)''{an)"' 
for all 771, n. It follows that 

Since the left-hand side does not depend on n and the right-hand side does 
not depend on m, this expression is a constant c independent of m and n: 

Pm = Poc™ and fj„ = o-qc". 

Since pm. and (T„ only have two possible values, we conclude that c = 1, and 
therefore that pm is independent of m and <t„ is independent of n. Therefore, 
for each x, 

7r(x + mui + nu2) = Ti{x)a'^ p^ , 

where p = X2 or 1, and cr = Ai or 1. (The choice of p and a can depend on x, 
but not on m or n.) □ 

We now use Proposition 2 to prove the main result of this section. 

Theorem 3. Assume that the translation invariant exclusion process 
with transition probabilities p{x,y) =p{y — x) has the property that there is 
no proper subgroup of Z'^ that contains P = {u ^ Z'^ -Piu) > 0}. If a is such 
that Ua is stationary, then there is a v & so that 7r(x) = 7r(0)e^^''"^ for all 
x. 
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Proof. Let ui,...,Uk be elements of P that span Z'^. To prove the 
theorem, it suffices to show that there is a G so that 

(20) {ui,v) =logpi, l<:i<k, 

where the quantities pi are the ones that appear in the statement of Propo- 
sition 2. By relabelhng, we may assume that ui,...,Ud is a basis for the 
vector space Q"^. Let v S R'^ be the unique solution of the equations (20) 
for 1 < i < d. For d <i < k, Ui can be written as a linear combination of 

Ul,...,Ud, 

d 

nui = ^ rijUj, 

where n and ni, . . . ,72^ are integers. Applying (19) twice gives 

7r(nnj) = 7r(0)/?" and 7r(niiii + • • • + n^n^) = 7r(0)p^^ . . . p^''. 
Therefore, 

Pi — Pi Pd ^ 

SO that 

d d 

nlog Pi = ^nj log pj = ^nj{uj,v) =n{ui,v), 

and (20) holds for this i as well. □ 

4. Product stationary measures on trees. For most of this section, we 
take S to be the homogeneous tree T in which every vertex has degree 3. 
(Similar examples with a general degree could be analyzed in a similar way, 
but we consider this case for concreteness.) Suppose q,r,s are positive, dis- 
tinct and sum to 1. We will use Theorem 1 to determine all stationary prod- 
uct measures for various exclusion processes on T that are homogeneous and 
nearest neighbor. We will initially consider p(x,y) so that the corresponding 
Markov chain has the property that, from each vertex, there are probabilities 
q, r, s of going out to the three neighbors, and also probabilities q, r, s coming 
in from the three neighbors. There are three inequivalent homogeneous ways 
of making these assignments that are given below as Cases 1, 2 and 3. The 
existence of nonreversible, inhomogeneous (i.e., with nonconstant density) 
stationary product measures for the process depends on the assignments. 

Case 1. The edge joining x and y is said to be of type 

1 if {p{x,y),p{y,x)} = {q,r}, 

2 a {p{x,y),p{y,x)} = {r,s}, 

3 if {p{x,y),p{y,x)} = {s,q}. 
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Assign labels 1,2,3 to the edges of T in such a way that for every vertex 
X the three edges incident to x consist of one of each of the three types. 
Up to isomorphism, there is only one way to do this. Once the labels are 
assigned, assign transition probabilities q, r, s in each direction for each edge 
in such a way that these assignments are consistent with the labels, and for 
each vertex, there are probabilities q,r,s of moving to the three neighbors. 
Again, up to isomorphism, there is only one way to do this. (In each case, the 
meaning of the word isomorphism should be clear. In the first occurrence, 
e.g., it means that, for any two assignments of labels to edges, there is a 
1-1 map from T with one assignment onto T with another assignment that 
respects the labels.) 

For a prescribed value of vr(xo) for a given vertex xq, there is only one 
choice of vr for which i/q, is reversible for the exclusion process with this 
kernel p{-); it is determined by 



for all neighbors x,y. We will show that the only stationary product mea- 
sures for this process are the homogeneous measures and the reversible mea- 
sures. 

Take any vertex x, and let xi,X2,X3 be the three neighbors of x, with 
the indexes chosen so that the edge joining x and Xi is of type i. If is 
stationary, then by Theorem 1(a), 

{q - r)l{^(^,^)=^(j,)} + {r- s)l{,r(x2)=7r(a;)} + {s - g')l{7r(x3)=7r(a:)} = 0. 

It follows that, for each x, either all or none of the neighbors Xi of x have 
the property that 7r{xi) = tt{x). However, it is clear that which of these 
situations occurs is independent of x, since any two neighbors must be in 
the same situation. Therefore, either vr is constant on T, or 7r(x) 7^ vr(y) for 
all neighbors x,y. In either case, vr is determined by its value at one site [by 
Theorem 1(b) in the latter case]. So, up to scaling, there are two stationary 
product measures, the first homogeneous and the second reversible. Note 
that these correspond to the two cases in Theorem 2.1 on page 380 in [13]. 

Case 2. The edge joining x and y is said to be of type 



Assign labels 1,2 to the edges of T in such a way that, for every vertex x, 
two of the three edges incident to x are of type 1 and the other is of type 2. 
Again, up to isomorphism, there is only one way to make the assignment. 
Assign transition probabilities to the edges in a manner consistent with 



TT{x)p{x,y) = TT{y)p{y,x) 



2 



1 



if {p{x,y),p{y,x)} = {q,r}, 
if {p{x,y),p{y,x)} = {s,s}. 



14 



M. BRAMSON AND T. M. LIGGETT 



the assignment of edge types. Each site x has exactly one neighbor y with 
p{x, y) = q, and one neighbor y with p{x, y) = r. 

To describe ah functions vr corresponding to stationary product measures, 
call a subset L of vertices of T a line if it is isomorphic to Z, and all edges 
joining vertices in L are of type 1. Then T can be partitioned into lines in 
a natural manner. Each edge of type 2 joins two lines. Since the transition 
probabilities are symmetric across edges of type 2, vr must take the same 
value at two vertices that are joined by an edge of type 2 [by Theorem 1(b)]. 
On each line, vr is either constant or of the form 



one sees this by applying Theorem 1(a), together with the isomorphism in 
the form L = {. . . , a;_i, xq, xi, . . .} with the appropriate orientation. So, the 
most general choice of vr so that Ua is stationary is obtained by (a) fixing 
the value of vr(x) at one vertex x, and (b) deciding for each line L whether vr 
is to be constant on L, or of the form (21). The stationary product measure 
is homogeneous if vr is constant on all lines, and is reversible if vr is of the 
form (21) on all lines. In all other cases, is neither homogeneous nor 
reversible. 

Case 3. The edge joining x and y is said to be of type 

1 if {p{x,y),p{y,x)} = {q,q}, 

2 ii{p{x,y),p{y,x)} = {r,r}, 

3 ii{p{x,y),p{y,x)} = {s,s}. 

Make consistent assignments of edge types and transition probabilities as 
before. In this case, the process is symmetric, so by Theorem 1(b), the only 
stationary product measures are the homogeneous ones, and they are, of 
course, reversible. 

In all three of the above cases, the homogeneous product measures are sta- 
tionary, but this is not necessarily the case for other homogeneous transition 
probabilities on the above tree T. For instance, one can choose the proba- 
bilities from a vertex to its three neighbors to be q, r, s as before, but with 
q = r ^ s (and, hence, s = 1 — 2q). Think of the tree as branching up, so that 
each vertex has two edges that go up and one that goes down. The probabil- 
ity q is assigned to each of the two upward edges and the probability s to the 
downward edge from each vertex. The corresponding Markov chain is not 
doubly stochastic. So by Theorem 1(a), there are no (nontrivial) homoge- 
neous stationary product measures, although the transition probabilities are 
themselves homogeneous. This, of course, does not occur for S = Z'^, where 
the Markov chain is doubly stochastic for such transition probabilities. 



(21) 
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The exponential growth of the number of neighbors within distance D of 
a given site in a tree can give rise to a phase transition for the corresponding 
stationary product measures. Consider, for example, the exclusion process 
in the preceding paragraph, but restricted to the (inhomogeneous) tree con- 
sisting of all vertices of T that are equal to or above a fixed vertex xq (with 
particles prevented from leaving this subtree). By Theorem 1, the stationary 
product measures for the process are reversible and satisfy 

7r(x) = {q/ (1 — 2g))"7r(xo) for D{xq,x) = n. 

When are such stationary product measures extremal stationary mea- 
sures? Theorem 2.1 in [7] gives, as a necessary and sufficient criterion [in the 
more general setting of countable S with irreducible that 

7r(x)/(l +7r(a;)) = oo. 

Here, this occurs for q/{l - 2q) G [1/2,2], that is, when q E [1/4,2/5]. This 
contrasts with S = Z, where none of the inhomogeneous stationary product 
measures is extremal, and with S = Z'^, d>2, where all stationary product 
measures are extremal. 

5. Stationary profile measures — necessary conditions for existence. In 

the next two sections we consider the exclusion process rjt on Z"^ with tran- 
sition probabilities p{x,y) = p{y — x) that satisfy 

< oo. 

X 

The analogues of blocking or profile measures in higher dimensions for such 
p(-) should be measures that are invariant under translations in certain di- 
rections, and have some specified limiting behavior in directions orthogonal 
to those directions. In this section we rule out the existence of stationary 
measures with these properties in certain directions. These results are ana- 
logues of the statement in one dimension that, if the particle drift is positive, 
and a stationary measure has density tending to zero in one direction and 
to one in the opposite direction, then it must be that the limit is one in the 
+00 direction and zero in the — oo direction. 

We begin with some definitions. Fix a nonzero vector v G R'^. li v ^ Z'^, 
we want to say that a measure /x is v-homogeneous if it is invariant under 
translations in all directions orthogonal to v. Unfortunately, this statement 
makes no sense for general v G R'^, since translations in directions orthogonal 
to V will, in general, lead outside of Z'^. In most cases, we will consider 
V £ Z'^, but to prepare for situations where v is more general (e.g.. Theorem 5 
below), we make a definition that is meaningful for all v € R'^, and agrees 
with the previous statement if v G Z'^. In general, a measure on {0, 1}^'' 
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will be called v-homogeneous if it has the following property: For each finite 
AcZ'^, there is a continuous function /a on R so that 

(22) ij{i]:i] = l onx + A} = fA{{x,v)), xeZ'^. 

It is easy to check that if t; G Z'^, this is equivalent to fi being invariant under 
shifts in all directions orthogonal to v, as we wanted. Note that the product 
measures Ua considered in Theorem 2 are i;-homogeneous for general v G R'^, 
since in that case, 

.Av:vix) = l} = ^^^-0:^. 

This supports the view that our definition is a good one in this context. Of 
course, a measure is u-homogeneous if and only if it is cu-homogeneous for 
any c£ R \ {0}, so there is no difference in statements of hypotheses between 
assuming v Z'^ and assuming v G Q'^. 

The continuity assumption on in the above definition is not an addi- 
tional requirement if u G Z'^, since then {x,v) only takes integer values. On 
the other hand, for general v G R'^, if we did not require /a to be contin- 
uous, property (22) would contain little, if any, information about The 
problem is that the map x — > {x,v) could easily be one-to-one, and in this 
case, any function of x is a function of {x, v). In remarks following the proofs 
of Theorems 4 and 6, we will indicate where the continuity property would 
arise if we were to take v G R'^ instead of v & Z'^. 

Note that if u G Z'^, for example, (22) does not determine the values of 
/a(s) for all s £ R. However, by choosing natural versions of /a, we may 
assume for general v G R'^ that 

(23) f^+A{s) = fA{s + {x,v)) foT xeZ^ and seR. 

An important property of i)-homogeneity is that the class of -y-homogeneous 
measures is preserved by the evolution of the exclusion process. This is clear 
if V G Z'^. To see it in general, let S{t) be the semigroup for the exclu- 
sion process, and for finite A C Z'^, let XAiv) be the indicator function of 
{ri:ri=l on A}. Then S{t)xAi'n) is continuous in so it can be uniformly 
approximated by functions of the form J27=i ^iXBi ■ Furthermore, since the 
exclusion process is translation invariant, S{t)xx+A is then automatically 
uniformly approximated by J22=i(^iXx+Bi, where the uniformity applies to 
both T] and x. Therefore, if /j, is u-homogeneous, 

nS{t){r]:i] = l on x + = J S{t)xx+Adfi 

is uniformly (in x) approximated by 

n „ n 

Xx+B,dfi = J2cifBi{{x,v)). 

i=l i=l 
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It follows that fiS{t){r} : r/ = 1 on x + A} is a continuous function of {x, v) as 
required. 

Unfortunately, the class of u-homogeneous measures is not preserved un- 
der weak convergence for general v because of the continuity assumption on 
/a- To see this, note that a product measure is u-homogeneous if and 
only if a{x) is a continuous function of {x,v). Suppose then that /„ is a se- 
quence of continuous functions on R that converges pointwise to a function / 
that has a jump discontinuity. Let an{x) = fn{{x,v)) and a{x) = f{{x,v)). 
If the coordinates of v are linearly independent over the rationals (so that 
{{x,v) :x ^ Z'^} is dense in R), then i/^n is t;-homogeneous, is not, yet 
converges weakly to z^q. As a consequence of this, we take v & Z'^ in. most of 
our results. 

A t^-homogeneous measure that is asymptotically equal to in the —v 
direction and asymptotically equal to 6i in the v direction, in the sense that 



will be called a v-profile measure. Our main objective in this section is to 
demonstrate the nonexistence of -y-profile stationary measures under suit- 
able assumptions on p{-) and v. In the next section we will provide some 
information about the structure of the set of -y-profile stationary measures 
in case of existence. 

The results of Section 3 give some examples of existence of ?;-profile sta- 
tionary measures. For instance, in the remark following Theorem 2, we con- 
sidered a general nearest neighbor process on Z^. Suppose for concreteness 
that pi > qi and P2> q2- Then, the four stationary product measures fi de- 
scribed there are, respectively: (i) homogeneous, (ii) v-profile with v = (1,0), 
(iii) f -profile with v = (0,1), and (iv) f -profile with 



This example illustrates why we do not want to restrict consideration to 
V ^ Z'^ (or, equivalently, to w G Q'^)- Even if the pi and qi are all rational, 
the above v need not be a multiple of a vector in Z'^. 

The proof of Theorem 4 below (and that of later results) uses the coupled 
process {'r]t,Ct) which has proved so useful in the study of the exclusion 
process. Each of the two coordinates is assumed to be a copy of the exclusion 
process, and the coupling is defined by saying that particles in the two 
processes move together as much as possible. This coupling has the property 
that sites at which rjt{x) ^ Cti^) can move and disappear, but they cannot 
be created. 

To describe the joint process in more detail, call a site x at which r/(x) 7^ 
({x) a discrepancy. There are two types of discrepancies, according to whether 



(24) 



lim /yi(s) = and lim /a(s) = 1 
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ri{x) = 1, <^(x) = or r]{x) = 0, C{x) = 1. Say that a site x with no discrep- 
ancy is of type 1 if r/(x) = ({x) = 1 and of type if 7]{x) = ((x) = 0. Consider 
now the possible transitions for the coupled process involving two sites x^y. 
If X and y are both of type 1 or both of type 0, nothing happens. If x is 
of type 1 and y is of type 0, the particles at x in both configurations move 
together to y at rate p{y — x). If there are discrepancies at both x and y, 
nothing happens if they are of the same type. If they are of opposite type, 
the particles in the two configurations move independently to the other site 
at the appropriate rates [one at rate p{y — x) and the other at rate p{x — y)]; 
after the transition, one of the sites x and y is of type 1 and the other is of 
type 0, that is, two discrepancies of the process have disappeared. Finally, 
if there is a discrepancy at x, and y is of type or 1, then the particle in 
the configuration for which the values at x and y differ moves to the other 
site at the appropriate rate. For more on this coupling and its application 
to the exclusion process, see Section 2 of Chapter VIII in [13]. 

The concepts ^-homogeneous and u-profile have natural analogues for the 
measures on {0, 1}^ x {0, 1}^ that arise as the measure of the coupled 
process. To define the first, the left-hand side of (22) is replaced by any 
cylinder probability, while for the second, one can simply say that on 
{0, 1}^ X {0, 1}^ is t;-profile if it is ^-homogeneous and each of its two 
marginals is f-profile. 

Our next result. Theorem 4, asserts that existence of f-profile stationary 
measures is only possible if v and the mean vector of p{-) have nonnega- 
tive inner product. Unfortunately, the proof requires that v be in Z'^. The 
problem that arises when v is general is closely connected to the continuity 
assumption that we made in the definition of ti-homogeneity. (For more on 
this point, see the remark following the proof.) 

The proof of Theorem 4 is based on the idea of the proof of Proposition 3.4 
in [1]. It involves a limiting argument that requires some domination. The 
following two lemmas provide the domination and main limiting statement 
that are needed. Neither involves the exclusion process or the kernel p{-). To 
begin, let vi, . . . ^ Z"^ be an orthogonal basis for R'^, with vi = v. Such 
a choice exists since the Gram-Schmidt orthogonalization procedure pro- 
duces vectors with rational coordinates if the original vectors have rational 
coordinates. For a positive integer n, let 

Bn = {xe Z'^:\{x,Vi)\ <n,l<i<d}. 

The statements of the two lemmas should be intuitively clear. (The precise 
constant in Lemma 1 is not important.) The reader should simply think of 
replacing the sums over points in Z'^ with integrals over R'^ in Lemma 1, for 
example. 
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Lemma 1. 

Y.\1bM - IbAu -z)\< 2(2n + 1)'^"^ Kz, v,)\. 

u j=l 

Proof. Write the left-hand side above as 
5]|lB„(n)-lB„(n-z)| 

u 

(25) 

The first term on the right-hand side is at most 
d 

(26) ^#{u:\{u,Vi)\ <n\fl<i<d,\{u,Vj) - {z,Vj) \ >n). 
i=i 

The mapping T from Z'^ to Z'^ defined by Tu = {{u,vi), . . . , {u,V(i)) is one- 
to-one (though not onto). Write u = {u^ ^),...,?x('^)) in Cartesian coordinates, 
and define 

C = {u:\u'^'^\ <nyi<i<d, \u^^^ - {z,Vj)\ > n} 

for fixed j and n. Then the jth summand in (26) is the number of elements 
in T~^C. This is at most the number of elements in C, which is, in turn, at 
most (2n -|- l)'^~^\{z,Vj)\. Summing on j and applying the same estimate to 
the second term on the right-hand side of (25), we see that (25) is bounded 
above by 

(27) 2(2n + l)'^-^X:K^'^^-)l' 
which completes the proof of the lemma. □ 

Lemma 2. There is a subsequence rij^ and a positive constant K so that 

(28) lim nl-'^i^{u G : {u, v)=r) = K 

fc— >oo 

for every r G F = {(u, u ) : n € Z'^} = aZ for some a > 0. 

Proof. We first show that n^~'^^{u G B^ '■ {u,v) = 0) is bounded away 
from and oo. As already observed, the mapping T defined in the proof of 
Lemma 1 is one-to-one, so 

#(u G Bn : {u,v) = 0) = #(u G Z'^:{Tu)^^^ = 0, |(rn)(^')| < n V2 < j < d) 

< G Z'^ : m(^) = 0, \u^^^\ < n V2 < j < d) 

= {2n + lf-\ 
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This gives the upper bound. For the lower bound, note that u = /c2f 2 + • • • 
+ kdVd S Bn ii ki^ Z and \ki\ < n|uj|~^ for 2 <i < d, so 

d 

#(n G Bn : {u,v) = 0) > n(2n|wr' " !)• 

1=2 

We may now choose the subsequence so that the hmit in (28) exists and 
is positive when r = 0. 

To complete the proof of (28), it now suffices to prove that, for r € T, 

lim ni-^[#(n G Bn : {u, v)=r)- #(u G Bn : {u, v) = 0)] = 

n — ^00 

(which is true along the full sequence). Since r G T, there is a u; G so that 
{w, v) = r. Therefore, 

#(u G Bn : (n, v)=r) = #{u e Bn'. {u - w,v) = 0) 
= #{u£Bn-w:{u,v) =0), 

from which it follows that 

#(n G Bn : {u, v)=r)- #(n G B„ : (n, v)=0)= 0{n''-^), 

as required, by the argument that led to (27). The only significant difference 
is that there is an additional constraint in the definition of the analogous 
C — the constraint that u^^^ = 0. This accounts for the reduction in the power 
of n from d — 1 to d — 2. □ 

Theorem 4. Take v G Z'^\{0} and suppose that fi is a v-profile measure 
that is stationary for the exclusion process. Then, (X^x ^P(^)' — 0- V Pi') 
is irreducible, the inequality is strict. 

Proof. Let u he a coupling measure that is f -homogeneous, is station- 
ary for the coupled process introduced above, and has first and second coor- 
dinates with measures /x and 1^1/2 j respectively. Such a measure exists. (See 
page 383 in [13], where this statement is proved without the part about v- 
homogeneity; the same proof applies to the f-homogeneous context because 
the property of u-homogeneity is preserved by the evolution of the exclusion 
process and by its coupled version.) The net rate at which the total number 
of discrepancies in Bn changes is zero, since the process is in equilibrium. 
On the other hand, transitions involving discrepancies of opposite type can 
only lower the number of discrepancies in Bn- Since /i is v-profile, there must 
be discrepancies of opposite type in Bn if n is sufficiently large. In the irre- 
ducible case, they will cancel each other out at a rate that is at least of order 
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^ Therefore, the net rate at which discrepancies enter i?„ is nonnegative 
in general, and is at least £n'^~^ in the irreducible case. Explicitly, 

^ p{x - y)u{{r], C) : r]{y) / C{y),v{x) = C{x) = 0} 

+ p{y-^y{{v,C)-v{y)7^({y),v{x) = C{x) = i} 

x£B„,y<^B„ 

(29) > pix-y>{iv,0-vi^)^Cix),viy) = Ciy) = ^ 

xeBn,yfBn 

+ piy - x)^{i^:C)-vix) j^cix), 

xeB„,y<^Bn 

r]{y) = C{y) = 0} + en^-\ 

where e > in general, and e > in the irreducible case. 

To take advantage of the ?;-homogeneity assumption, define functions gz 
and hz by 

9z(.{x, v)) = iy{{ri, C) : r/(x + z) C{x + z),v{x) = C{x) = 0}, 
hz{{x,v)) = u{{ri, C) -Vix + z) C{x + z),v{x) = C{x) = 1}. 

Then, (29) becomes 

Y P{x-y)gy-x{{x,v))+ Y P{y-x)hy^^{{x,v)) 

x&B„.y(^Bn xeBn,yfBn 

> Y pix-y)hx-yi{y,v)) 

xeBri,yfBn 

+ Y p{y-x)g^-y{{y,v)) + en'^-^. 
x€Bn,y<^B„ 

Making the changes of variables z = x — y or z = y — x and u = x or u = y 
in the four sums, and noting that, for any sets B and C, 

this leads to 

Ypiz)Y9-zi{u,v))[lBM - ^B„{U-Z)] 
z u 

(30) 

> Ypi^) Y ^-(("' ^)) [l^n (u + z)- 1b„ (n)] + en'^-i . 

z u 

Next, we will divide (30) by n'^~^ and pass to the limit (along a subse- 
quence) using the dominated convergence theorem. The domination is pro- 
vided by Lemma 1. Using this and Lemma 2, we will prove that, for any 
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function g on F that satisfies 

lim g{l) = 1 and lim g{l) = 0, 

I—*— CO /— ►oo 

(31) lim nl-''Y.g{{u,v))[lB^^{u)-lB„Ju-z)]=K{z,v)/a. 

u 

Any such function g can be uniformly approximated by functions of the form 

m 

where each Ij G T and Sj=i = 1. By Lemma 1, it suffices to prove (31) for 
a function g of the form g{l) = IjKr}) where r G F. In this case, if {z, v) > 0, 
the left-hand side of (31) becomes 

(32) lim 4~^#(n eB^^-.r- {z, v) < {u, v)<r), 

fc— >oo 

so that (31) follows from (28), since there are {z,v)/a elements of T in the 
interval (r — {z,v),r]. A similar argument gives (31) when {z,v) < 0. 
Note that by (24), 

lim g^{l) = i, lim g^{l) =0, 

t— oo I— too 

lim ft,2(/) =0, lim h^{l) = ^ 

Z— >— oo l^oo 

for every z. Therefore, we may apply (31) to g{l) = 4:gz{l) for any z. The 
conclusion is that the limit of the left-hand side of (30) divided by n'^~^, 
along the sequence n^, is ^J2zPi^){^^'^) ■ Similarly, applying (31) to g{l) = 
1 — 4:hz{l), it follows that the limit of the right-hand side of (30) divided by 
n'^~^, along the sequence Uf^, is e — ^J2zP(^)i^^'^)- Therefore, (30) implies 
that J2zPi^){^^''^) — '^cle/K, as required. □ 

Remark. As mentioned earlier, the main difficulty in extending The- 
orem 4 to general v is related to the continuity assumption on Ja in the 
definition of ?;-homogeneity. In the proof, we used a w-homogeneous cou- 
pling measure with prescribed f-homogeneous marginals. If v were general, 
the continuity assumption would be needed in the treatment of (31) where 
we approximated g by step functions. The issue is the existence of a coupling 
measure that satisfies the continuity property. In the proofs of analogous re- 
sults in [13] (see pages 143 and 383), one starts the coupled process with 
the measure that is the product of the two given marginals, and passes to 
a Cesaro limit of the measure as t — > oo. What is to guarantee that this 
limiting measure satisfies the needed continuity property? Perhaps a sta- 
tionary measure for the marginal or coupled process automatically satisfies 
u-homogeneity with the continuity property. We do not know whether this 
is always the case. 
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Under a partial symmetry assumption, it is not difficult to prove a similar 
result that allows v G R'^: 

Theorem 5. Take v G and suppose p{z) =p{—z) for all z for which 
{z,v) 7^ 0. If n is V -homogeneous and stationary for the exclusion process, 
then ii{rj : rj{x) = 1} is constant on the smallest subgroup of Z'^ containing 
P = {x € Z'^:p{x) > 0}. Therefore, there are no v-profile stationary mea- 
sures for the process in the irreducible case. 

Examples to which this result can be applied are the nearest neighbor 
exclusion processes in which the transition probabilities are symmetric in 
all directions but one, say ei, and with a drift in the direction of ei. The 
conclusion in these examples is that there cannot be any Cfc-homogeneous 
stationary measure /i with nonconstant density if fc > 1. 

Proof of Theorem 5. Applying the generator of the exclusion process 
to the function rj ^ vi^): integrating with respect to fi, and using station- 
arity, gives, for each x ^ Z'^, 

^p{y - x)n{ri : ri{x) = 1, r?(y) = 0} 
y 

(33) 

= ^p{x- y)n{r] : r]{x) = 0, ry(y) = 1}. 

y 

Making the change of variables z = y — x in the sums and using (22) yields 

Y^p{z)[f{{x,v))-f.{{x,v))] 

(34) 

z 

where we have used the shorthand 

/(s) = f{o}{s) and /^(s) = /{o,4(s)- 
Rewriting (34) gives 

(35) - /^(^)] = Ep(^)[/(« - ")) - /^(« - ^))] 

z z 

for all s G r = {(x, f ), X E Z'^}. In obtaining (35), we have used the relation 

(36) /_^(s) = /{o,-4(s) = f~z+{o,z}{^) = f{o,z}{^ ~ '^)) = ■/'^(^ ~ i^^v)), 

where we have used (23) in the third equality. So far, we have used prac- 
tically none of the assumptions of the theorem. We point this out because 
(35) and (36) will be used later in the proof of Theorem 7. 
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Now, by the symmetry assumption on p{-), the terms in (35) involving 
cancel out, since 

Y^p{z)[f,(s-{z,v))-Us)]= Yl Piz)[f^.{s) - Ms)] 
(37) = J2 Pi-z)[f-z{s)-fz{s)] 

z: (z,t>>^0 

= pmz{s) - f-z{s)]- 

z: {z,v)j^O 

We have used (36) in the first step, the symmetry assumption in the sec- 
ond and replacement of z by —z in the third. But the second and fourth 
expressions in (37) are negatives of one another, so they must be zero. 
Using the fact that the sums in (37) vanish, (35) leads to 

fis) = Ypiz)f{s-{z,v)), sev. 

z 

This says that g{x) = f{{x,v)) is a bounded harmonic function for the ran- 
dom walk on Z"^ with probabilities p{—z) of moving from x to x + z. It 
follows that / is a constant on the irreducible classes for this random walk 
by the Choquet-Deny theorem — see [3], for example. □ 

6. Properties of stationary profile measures. In this section we study 
exclusion processes on Z'^, d>2, whose transition probabilities are transla- 
tion invariant and have a finite mean. Assume, in addition, that the process 
is irreducible. We begin by using techniques similar to those in Sections 
2 and 3 of Chapter VIII in [13] to determine the structure of f-profile sta- 
tionary measures when they exist. Throughout this section, v will be a fixed 
nonzero element of Z'^. 

Theorem 6 below has two statements. We will prove the first in detail, 
but will only sketch the proof of the second part. Three applications of 
Theorem 6 follow its proof. In the statement below, the inequality /xi < 112 
refers to stochastic monotonicity. See page 71 in [13] for the definition. Recall 
that an equivalent statement is that fii and H2 can be coupled so that any 
site at which the fii configuration has a 1 must also have a 1 in the /U2 
configuration. 

Theorem 6. Suppose that fii and ^2 o,re extremal stationary measures 
for the exclusion process. 

(a) // III and fi2 olt^ v-profile, then either p,i < fi2 or ^2^ P-i- 

(b) If Hi and fi2 are v-homogeneous and {J2z zp{z),v) =0, then either 
fj.1 < H2 or ^2 < 
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Proof, (a) It is enough to prove that every extremal stationary v-profile 
measure ly for the coupled process concentrates on configurations (r/, (^) that 
have no discrepancies of opposite type, that is, that each such i/ satisfies 

C) : vix) = C(y) / viy) = Cix)} = 

for all X, y G Z"^. To do so, define -Bm,n for 1 <m <n slightly more generally 
than the corresponding sets Bn used in the proof of Theorem 4: 

Bm,n = {x G Z'^ : \ {x,vi) \ < m, \ {x,Vi) \ < n for 2<i < d}, 

where vi,...,V(i G are orthogonal and vi = v. [Later in the proof, we 
will take m,n — > oo with m = o{n).] We will consider the contributions to 
the expected (with respect to i^) net rate at which discrepancies in i?m,n 
disappear. The expected total net rate is zero because the coupled process 
is in equilibrium. 

If discrepancies of opposite types appear with positive u probability, then 
they will disappear in Bm,n at a rate that is at least of order n'^~^, since (by 
f-homogeneity) there will be at least this many pairs of sites x, y in -Bm,n 
with the same difference x — y and the same (by irreducibility) nonzero 
probability of having a discrepancy of one type at x and of the other type 
at y. Therefore, it suffices to show that the net rate at which discrepancies 
enter Bm,n is o{n'^~^). 

The idea is that discrepancies can enter i?m,n across two boundary faces 
on which the density of particles is close to in one case and close to 1 in the 
other (by the 7;-profile assumption), or across the other boundary faces. In 
the first case, there are roughly n'^~^ locations to consider, and each makes 
a contribution to the rate of entry of discrepancies that is o(l), so the total 
contribution is o{n'^~^). In the second case, the number of locations is of 
order mn'^"'^, which is o(n'^~^) if we choose m = o{n). These correspond to 
cases (B) and (A) below, respectively. 

To carry out the details, consider sites x G Bm,n,y ^ Bm,n, and the ex- 
pected rate at which a discrepancy moves from x to y or from y to x. For 
example, the expected rate at which a discrepancy moves from x to y is 

p{y-x)u{{7],c)-v{x) 7^C{x),v{y) = C{y) = 0} 

+ p{x-y)v{{ri,Q:r]{x) ^ (;{x),r]{y) = ({y) = 1}, 

with a similar expression for the expected rate at which a discrepancy moves 
from y to x. Sums of these expected rates over appropriate choices of x and y 
will be bounded in the various cases. 

Since y ^ Bm,n, either \{y,vi)\ > m or \{y,Vi)\ > n for some 2 < i < d. 
There are two cases to consider: 



(A) \{y,Vi) \ > n for some i> 1. 
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(B) \{y,vi) \ > m and \ {y,Vi) \ < n for alH / 1. 

The contributions to the expected rate of motion of discrepancies involv- 
ing all x,y in case (A) (corresponding to motion from points x inside the 
box to y outside the box or vice versa) for a given i are bounded above by 
(letting z = y — x) 

Y^[p{z) +p{-z)]#{x:\{x,vi)\ < m,\{x,Vj)\ <n\/j> l,\{x + z,Vi)\ > n). 

z 

This sum is at most 

Cmn'^~^^\z\p{z) 

z 

for some constant C. This is o(n^~^), provided we let m,n — > oo with m = 
o(n). 

Next, we consider the total contribution to the expected rate at which 
discrepancies enter Bm,n for x and y in case (B). Letting z = y — x, we see 
that this total expected rate is bounded above by 

^[p{z)+p{-z)] u{{r],C)-vix + z)^C{x + z)}. 

\(x,Vi) + {z,Vl)\>m 

The number of terms in the inner sum is bounded by a constant multiple of 
|^;|n'^~^, uniformly in n,m and z. Therefore, by the dominated convergence 
theorem and the first moment assumption on p(-), the above expression will 
be o(n'^~^), provided that, for each z, the ratio of the inner sum to the 
number of summands in the inner sum tends to zero as m, n — > oo. Consider 
the summands for which {x,vi) + {z,vi) > m, for example — those for which 
{x, vi) + {z, vi) < —m are handled in a similar manner. For these summands, 
we have m— {z,vi) < {x,vi) <m, and 

^{{Vj C) : v{x + ^) / C{x + z)} < ni{ri : ri{x + z) =0} + fi2{v ■ vi^ + = 0} 

= 2 - f{z},l{{x,v)) - f{z},2i{x,v)), 

where we have used fA,i for i = 1 and 2 to denote the function defined in (22) 
corresponding to the measure //j. Since fii and /X2 are u-profile measures, the 
right-hand side above tends to zero as m — > oo. (Recall that v = vi.) This 
completes the consideration of case (B). 

Sketch of proof of (b) . The only part of the above argument that uses the 
assumption that the measures Hi are u-profile, as opposed to f -homogeneous, 
is the treatment of case (B). In that context, we used the fact that the density 
of discrepancies near the faces of Bm,n that are orthogonal to v is small when 
m is large because the density of O's (for the face in the +v direction) or I's 
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(for the face in the —v direction) is small. The point is that, at a discrepancy, 
one coordinate must be and the other must be 1. If we only assume that the 
measures are f -homogeneous, the net expected rate at which discrepancies 
cross the faces of Bm,n orthogonal to v has to be shown to be small for a 
different reason. This reason is, of course, related to the assumption that the 
mean vector for p{-) is orthogonal to v. The argument parallels the proof of 
Theorem 3.14 on page 391 in [13]. 

Here is the main idea. In a Cesaro sense (with respect to shifts in the zizv 
directions), the coupling measure v, near the faces of Bm,n that are orthog- 
onal to V, is nearly translation invariant. It is translation invariant in the 
directions orthogonal to v by the assumption of f-homogeneity, and in the v 
direction as a result of the Cesaro averaging. So, in the limit as oo, the 
measure is both stationary for the coupled process and translation invari- 
ant. This implies that in this limit, the marginals fxi and //2 are stationary 
for the exclusion process and translation invariant. But this implies that they 
are exchangeable by Theorem 3.9 on page 388 in [13]. Furthermore, they are 
coupled together so that there are no discrepancies of opposite type. 

By decomposing the limiting coupling measure according to the type of 
discrepancy that occurs, we find that we are essentially in the following 
situation. We can assume that only one type of discrepancy occurs, so that 
the measures //i and /i2 are exchangeable, and are coupled by a measure v 
that satisfies 

K(r/,C):??(x) = l,C(^)=0}=0 

for all X. We will use a suggestive notation for probabilities related to i'. For 
example, we will write i^{(r/, () : rj{x) = 0, rj{y) = C,{x) = ({y) = 1} as 

.0 1, 

Entries on the left correspond to site x, while entries on the bottom corre- 
spond to configuration rj. The total expected net rate at which discrepancies 
of this type enter Bm,n across the face orthogonal to v in the +v direction 
in case (B) is then 



(y,vi)>m 
\(y,vj}\<n VjVI 



+ p{y 



1\ /II 
1 



1 1\ /I 

1 -'^ 



Using the fact that < a.s. with respect to z/, the differences that appear 
above can be written as 



1\ /II 
/ 



(o l) =MC:C(x)=0,C(y) = l}-/^i{r?:r?(x)=0,r?(y)=l} 
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and 

^(i o)"^(o o)=MiO?:'?W = i.^(y) = o}-MC:C(^) = i,C(?/)=o}. 

By exchangeability of the marginals, the right-hand sides above are inde- 
pendent of X and y and are negatives of one another, so we will call them 
A and —A, respectively. It follows that the total expected net rate at which 
discrepancies enter Bm,n, across the face orthogonal to v in the +v direction 
is, in case (B), 

A [p{x - y) - p{y - x)]. 

x£Bm,n,{y,vi)>m, 
\{y,Vj)\<n Vjyi 

Except for the factor of A, this can be written as 
Yp{z)[#{x £ Bm^n ■■ {x,vi) > m + {z, vi) ,\{x , Vj) - {z,Vj)\ <nyj^ 1) 

z 

- #(x G Sm,„, : > m - (z, vi), fj) (z,Vj)\ < n Vj / 1)]. 

It is not hard to check that this expression is asymptotic to a negative 
multiple of 

n'^'^Yp{z)(,z,vx). 

z 

Since the mean vector of p(-) is orthogonal to wi, this is o(n'^^^), as required. 
□ 

Remark. If one tried to prove Theorem 6 for general v G i?*^, one would 
need the continuity assumption in the definition of u-homogeneity for the 
coupling measure in order to know that, if v concentrates on configura- 
tions with discrepancies of both types, then the rate of loss of discrepancies 
inside -Bm,n is of order . As explained earlier, we do not know whether 
the continuity assumption on the marginals [li is inherited by the coupling 
measure v. 

As a consequence of Theorem 6, we can explicitly identify all stationary 
ti-profile measures when (6) is satisfied. 

Corollary 1. Suppose that (6) holds, and let 
ac[x) = -. — r, c> 0. 

Then, the set of all extremal stationary v-profile measures is exactly {v^c j c > 
0}. 
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Proof. Let // be any extremal stationary f -profile measure. The mea- 
sure Voic has the same properties: it is stationary by Theorem 2 and it is 
extremal by Theorem 2.1 of Jung [7], since 

^ adx) [1 - ac{x)] = oo. 

X 

(Jung's context assumes reversibility, but this property is not needed in his 
proof.) Therefore, by Theorem 6, for each c > 0, either /_f < i^a^ or z/q^ < fj,. 
Since the measures z^^^ are stochastically monotone in c, there is some < 
c* < oo so that /i < z^oc for c > c* and fi > for c < c* . Since the family 
i^a^ is weakly continuous in c, it follows that ji = i^a^* ■ D 

Corollary 2. Suppose that 

xp{x),v^ = 0. 

Then, the extremal stationary v -homogeneous measures are exactly the ho- 
mogeneous product measures. 

Proof. The proof is the same as that of Corollary 1. The role that 
played in the proof of Corollary 1 is now played by the homogeneous product 
measures. In this case, Theorem 1.17 on page 216 in [14] can be used as an 
alternative to Jung's theorem. □ 

Finally, we will combine Theorem 6 with part of the proof of Theorem 5 
to show that, under an additional assumption on p(-), all -y-profile stationary 
measures have a property analogous to that of a blocking measure in the 
sense of Bramson, Liggett and Mountford [1]. The proof is based on the 
proof of Lemma 6.4 of that paper. The simplest way to explain the idea of 
the proof of Theorem 7 is to refer to identity (6.12) in the statement of that 
lemma. It states that if d = 1 and z^ is a stationary measure for the exclusion 
process with transition probabilities p{-) that has a finite expected number 
of ones to the left of the origin and a finite expected number of zeros to the 
right of the origin, then 

oo oo 

n^p^—n) = ^ nM(n)[p{n) — p(— n)], 

n=l n=l 

where M{n) is the i^-expected number of sites x at which there is a one at 
X and a zero at x -\- n. This provides an a priori bound for the quantitites 
M{n) in terms of the second moments of the negative jump probabilities 
for the particles. The proof below develops this idea for systems in higher 
dimensions. Most of the work is aimed at proving the convergence in (39). 
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Theorem 7. Suppose that 

(a) p{z) >p{—z) whenever {z,v) > 0, 

(b) p{z) > p{—z) for some z such that {z,v) > 0, 
and 

(c) Ez: (z, v)>0 {Z, V)^ <00. 

If fi is a V -profile extremal stationary measure for the corresponding exclu- 
sion process, then fi is concentrated on configurations rj that satisfy r]{kv) = 
for all sufficiently large negative k and i]{kv) = 1 for all sufficiently large 
positive k. 

Proof. We begin by restating (35), which only requires the assumption 
of stationarity. In the present case, F = NZ for some integer > 1, since 
V ^ Z'^. We therefore have 

Y,p{z)[f{kN)-f,{kN)]=Y,p{z)[f{kN-{z,v))-f,{kN-{z,v))]. 

z z 

Sum this identity over m < k < n. For fixed z, make the change of variables 

....<^ 

in the sum over k of the terms on the right-hand above. A significant amount 
of cancellation occurs with terms on the left-hand side, whose precise nature 
depends on the sign of {z,v). The result is 

n 

Y: P{z) Yl [fikN)-f,ikN)] 

z:{z,v)>0 k=n+l-({z,v)/N) 
n-{{z,v)/N) 

- E P(^) E [f{kN)-f,{kN)] 

z:{z,v)<0 k=n+l 

m 

= E P(^) E [f{kN)-f,{kN)] 

z:{z,v)>0 k=m+l-{{z,v)/N) 
m-{{z,v)/N) 

- E Pi^) E [f{kN)-fzikN)]. 

z : {z,v)<0 k=m+l 

Since the left-hand side above becomes the right-hand side if n is replaced 
by m, it follows that the left-hand side is independent of n. Since /i is v- 
profile, this expression tends to zero as n — > ±00. Therefore, it is zero for all 
n. In other words, 

n 

Y P{z) E [f{kN)-f,{kN)] 

z:{z,v)>Q k=n+l-{{z,v)/N) 
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n~({z,v)/N) 

E Pi^) E [fikN)-UkN)] 

z:{z,v)<0 k=n+l 



over all n. Next, sum this identity for —M <n < M. The result is 

k '-^ 



E p{z)Y,[f{kN)-UkN)] 

z: {z,v)>0 



: E p{z)Y.[fikN)-UkN)] 

z : (z,v)<0 k 



^ AM-(-M) V/c 
{z,v) 



N 

kAM - {k + 



N 



V (-M) 



Replace z by —z on the right-hand above, then make a change of variables 
k ^ k + {{z,v)/N), and finally use (36). The right-hand side above then 
becomes 



E p{-z)Y.[f{kN+{z,v))-MkN)] 

z:(z,v)>0 k 



k + 



N 



^ AM-{-M)\Jk 



Therefore, moving some terms from the resulting right-hand side to the 
left-hand side, we obtain 



(38) 



where 



E \p{z)-p{-z)]Y.yikN)-f^(kN)]^M{k,z) 

z : {z,v)>0 k 

= E p{-z)Y.[fikN + {z,v))-f{kN)]cM{k,z), 

z:{z,v)>0 k 



CM{k,z) 



(^k+^^^ AM-{-M)Vk 



By comparing fi with a translate of fi in the v direction, we see from 
Theorem 6 that f{kN) is a monotone function of k. Since /x is a f-profile 
measure, it must be an increasing function of k. Therefore, all the summands 
on the right-hand side of (38) are nonnegative. Since 



CM{k,z) < 



N 



for all choices of the arguments, the right-hand side of (38) is at most 



N 



-1 



E p{-z){z,v)J2[f{kN + {z,v))-f{kN)] 



z:{z,v)>0 k 
2 : {z,v)>0 
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which is finite by assumption. In the identity above, we have used the fact 
that {z,v) is a multiple of N . 

The summands on the left-hand side of (38) are also nonnegative by as- 
sumption. Since 



lim CM^k, z) = 




we may apply Fatou's lemma to the left-hand side of (38) to conclude that 

z : {z,v)>0 k 

By assumption, there exists a z with {z,v) > so that p{z) >p{—z). For 
this z, it follows from this display that 

(39) Y.if{kN)-f,{kN)]<<^. 

k 

We claim that the convergence of the series in (39) does not depend on z, 
and hence it will converge for all z. Note first that by the definitions of 
/ and fz, 

(40) /i{r/ : r/(x) = 1, ri{x + z) = 0} = /((x, v)) - f,{{x, v)). 

To see that the convergence of (39) holds for all z, it suffices by (40) to check 
that, for any z and w, there exists an e > so that, for all x, 

fi{r] : r]{x) = l,r]{x + z) = 0} > efi{7] : r]{x) = l,r]{x + w) =0}. 

For this, it is enough to show that such an e exists so that 

li{ri : r]{x) = 1, r]{x + z) = 0, r]{x + w) = 1} 

> efj,{r] : r]{x) = 1, r]{x + w) =0, rj{x -|- z) = 1} 

for all X. To check this, since /i is stationary, it is enough to use irreducibility 
to find a path zq, zi, . . . , z^ with zq = z, z„ = w, and p(zj+i — Zj) > for 
each i, and then observe that there is a way to go from any configuration rj 
such that r]{x -|- z) = 1, r]{x + im) = to r]x+z,x+w (the configuration obtained 
from ?7 by interchanging the states at x + z and x + w) by successively moving 
particles along the path without interference from sites not on the path 
(shifted by x), so that the probability is positive that these transitions occur 
in the given order by time one, and that by time one, no other transitions 
involving sites on this path are attempted. Schematically, suppose rj has the 
following form along the shifted path: 

110 110 

and we need to move to 

10 110 1. 
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The way to do this is to move the rightmost 1 two steps to the right, the 
next 1 one step to the right, the next 1 two steps to the right, and finally 
the leftmost 1 one step to the right. 

We have shown that the series in (39) converges for all z and, in particular, 
for z = V. Recalling that |f p is a multiple of N, say, |wp = IN, we conclude 
that 

Y,fi{rj:v{kv) = + l)v) = 0} = $][/(MA^) - MklN)] < oo. 

k k 

Therefore, 

lim rj(kv) and lim rj{kv) 

fc— >— oo fc— ^+00 

exist a.s. Since // is a f-profile measure, these limits must be and 1, re- 
spectively, as claimed. □ 

7. Convergence to stationary profile measures. In previous sections we 
studied the stationary measures of exclusion processes. Here and in the next 
two sections, we give convergence results for exclusion processes with non- 
stationary initial measures. The exclusion processes are assumed to be on 
Z'^, d>2, with transition probabilities that are translation invariant, have 
finite mean and are irreducible. The main result in this section is Theorem 8, 
where we show that, given a one-parameter family of extremal stationary 
i;-profile measures, an exclusion process with a u-profile initial measure con- 
verges weakly to an average of these stationary measures, provided one is 
given certain bounds on the tails of the measures. As in Section 6, we assume 
u is a nonzero element of Z'^. 

In order to show Theorem 8, we will employ four lemmas. For the first of 
these, we introduce the following notation. Let vi,...,Vfi be an orthogonal 
basis for R'^ with Vi G Z"^ and vi =v; as mentioned before Lemma 1, such a 
basis can be constructed using the Gram-Schmidt procedure. We set 

L{z) = {x e Z'^ -.0 < {x - z,Vi) < {vi,Vi), 2<i<d} 

and 

L- (z) = L{z) n {x £ Z'^ : {x - z,vi) < 0}, 

L+{z) = L{z) n{x£ Z'^:{x- z, vi) > 0}. 

That is, L{z) is an infinite strip in the vi direction which has width \vi\ in 
the perpendicular directions, and L^{z) and L^{z) are the corresponding 
semi-infinite strips. For k = {k^'^\ . . . , k^'^'>) G Z^-i, set 

Lk{z)=Ll^z + Yj,^^v)j. 
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It is easy to see for given z that Lk{z), k £ Z'^ ^, partitions Z'^ and 

d 



(41) xeL(z), yGLkiz) 



d 

X 



i=2 
d 



y-Y^k^\,eL{z). 



i=2 

For a f -homogeneous measure A, we also set 

q{z) = E\#x(^L-{z):ii{x)=ll 

(,+ {z) = E^[#x G L+{z):r]{x) = 0]. 

In the following lemmas and in Theorem 8, we will employ f -homogeneous 
measures A satisfying 

(42a) q{0)<oo 

and 

(42b) £^{0)<oo. 

Note that such u-homogeneous measures are v-profile. 

We will repeatedly use the coupling = {i]t,Ct) introduced in Section 5 
for the exclusion process. Unless stated otherwise, we will assume that the 
initial measures of the two coordinates are independent. Let ji and A denote 
the initial measures of the coordinates, respectively, and 7 the corresponding 
product measure. 

Let be such a coupling, where fj, and A are t^-homogeneous measures 
satisfying (42b) for n and (42a) for A. The first part of Lemma 3 states 
that the expected number of discrepancies of type (0, 1), along a strip L{0), 
is nonincreasing in time. The second part of the lemma uses this to obtain 
uniform upper bounds in time on the tail of {z) for a measure A satisfying 
(42a), when a f -homogeneous stationary measure satisfying both parts of 
(42) exists. 



Lemma 3. Assume that /x and A are v -homogeneous measures satisfying 
(42b) and (42a), respectively. Then, 

(43) E'^y^x G -^^(0) : rit{x) = 0, Ct{x) = 1] is nonincreasing in t. 
Assume, moreover, that fj, is stationary and satisfies (42a). Then, 

(44) £l^{z)^0 as {z,v)^ -co 
uniformly in t. 
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Proof. Discrepancies of type (0, 1) can disappear, but cannot be cre- 
ated. So, in order to show (43), it suffices to sliow tliat tlie expected rate at 
wliich discrepancies leave L{0) equals the expected rate at which they enter 
at a given time, if one ignores terms corresponding to their disappearance. 
Let r{x,y;(^) denote the rate at which a discrepancy of type (0,1) moves 
from X to y; this depends on just p{y — x) and the present state ^. The 
expected rate at which such discrepancies leave L{0) at time t is 



^7 



x&L{0)y^L{0) 



fc^O 



Y Y ^i^^y^^t) 

xeL(0)yeLfc(0) 



which is finite because of assumptions (42b) and (42a) on fi and A. Since /x 
and A are f -homogeneous, so is 7, and hence so is the distribution of ^f. It 
therefore follows from (41) that the above sum equals 



YE- 



Y Y ^(^'2/;^* 

x-GL_fc(0)-yGL(0) 



Y Y ^i^^y^^t) 

Lx^L(0)yeL(0) 



which is the expected rate at which these discrepancies enter L{0). 

We now demonstrate (44). Let „;U denote the shift of /i by —nv and „7 
the product of „/U and A; since fi is ?;-homogeneous and satisfies (42b), the 
same is true for It follows from this and (43) applied to „;U that, for 
given e > and sufficiently large n, 

E'^^[#xeL{0):rit{x)=0, Ct{x) = 1] < e 

uniformly in t. 

On the other hand, since /i is stationary and is assumed to satisfy (42a), 

E"''[#x G L-{-2nv):r]t{x) = 1] = E^[#x G L'i-iiv) :7]t{x) = 1] < e 

for sufficiently large n and all t. Together with the previous inequality, this 
implies 

q^{-2nv) = E^li^x £ L-{-2nv):Ct{x) = 1] < 2e 

for all t, since L~{—2nv) C L{0). This implies (44) for z = —2nv. The limit 
also holds for general z since L~{z) is contained in a finite union of appropri- 
ate L^,{—2nv), j = 1,... ,2'^"^, where the individual coordinates of kj differ 
by at most 1 and n is chosen so that 

-2{n + l){v,v) < {z,v) < -2n{v,v), 

and since for any choice of k, 

q(-2nv + Yk^'^Vi] =^A,(-2™)- 



i=2 



□ 



36 



M. BRAMSON AND T. M. LIGGETT 



Let 

L^(0) = L(0) n{xeZ'^: -N{v, v) < {x, v) < N{v, v)}, 

where N € Z+ , and let L^(0) = L(0) \ L^(0). Also, let 

BN = {xe Z'^:-N{vi,Vi) < {x,Vi) < N{vi,Vi) for i = 1, . . . , d}, 

N G Z~^, where {vi, . . . ,Vd} is the orthogonal basis for R!^ defined earlier, 
with vi = V. Lemma 4 provides upper bounds on the number of discrepan- 
cies for the coupled process = {r]t,Ct), starting from u-profile measures ji 
and A satisfying (42), where ^ is stationary. The first part gives bounds on 
discrepancies of types (0,1) and (1,0) existing simultaneously in Bn, and 
the second part gives bounds on the total number of discrepancies in L^(0). 

Lemma 4. Assume that fi and A are v-profile measures satisfying (42), 
and that fi is stationary. Then, 

(45) {discrepancies of opposite types exist for {ijtXt) inBj\f}^0 
as i —> oo for each N . Also, 

(46) E'^[# discrepancies of {r]t, (t) in ^^(0)] — > as N ^ oo 
uniformly in t. 

Proof. In (45), we may assume that is large. Let L^^b denote the 
union of the (4A^)'^~^ translations of L{0) containing points in B2n- It suf- 
fices to show for each t at which the probability in (45) is at least e > 0, 

ii^"^ [^discrepancies of in L7v,_b] — -E''''[#discrepancies of ^t+i in -^Af.s] < 

where e' > depends on e and N, since this implies there are only a finite 
number of such t at least distance 1 apart. 

As in the demonstration of (43), one can show that the expected number 
of discrepancies in L^^b is nonincreasing at each t by calculating the terms 
corresponding to their movement, but ignoring those corresponding to their 
disappearance. It therefore suffices to show that, for any configuration with 
a pair of discrepancies of opposite types at sites z,w € B]\f, there is a uniform 
positive lower bound on the probability of these discrepancies meeting each 
other or other discrepancies of opposite type after an additional unit of time, 
while remaining in B2N C L^^b- This is shown in Lemma 3.1 in [1]. The 
argument uses the irreducibility of p(-), and is similar to that near the end 
of the proof of Theorem 7, where one constructs a path zq, 2:1, . . . , z„, with 
zq = z, Zn = w, and p{zi^i — Zi) > for each i, along which the discrepancies 
may move while remaining in B2n, until meeting one another or some other 
discrepancy of opposite type. 
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To show (46), we apply (44) both to A and to the corresponding measure 
with the roles of and 1 reversed, to get 

£^^(— z)— >0 and £^^(z)— >0 as(2;,f)^oo 

uniformly in t. Since fi is stationary and is assumed to satisfy (42), fit satisfies 
the analogous limits. Together, these four limits imply (46). □ 

The following elementary lemma compares extremal stationary measures 
for the exclusion process on the sets 

DM={{r],0-Viz)>C{z) for zGSm}, 

where M G Z+. 

Lemma 5. Let 7 be any measure with marginals ^u^^ and fj,°'^ , where /x"^ 
and fi"^ are extremal stationary measures with /i^^jry : ry(0) = 1} < fj,"^{rj: 
r/(0) = 1}. Then, for any e > 0, 

(47) 7(^A/) < e 

for large enough M not depending on the choice of ^. 

Proof. Assume, on the contrary, that there exist such measures 7*^ and 
sets DM^, with ^^{Dm^,) > £ and — > 00, as fc — > 00. Choose a subsequence 
kn along which 7^^" ^ 7°° as n ^ 00 for some measure 7°° ; then 7°° has the 
same marginals Let 

D = {(7?, C) : r]{z) > C{z) for all z}; 

we claim that 7°°(D) = 0. It will follow that 7^" (Da/) < e for large enough 
n and M, which will contradict our assumption on 7''. 

To see that ^°°{D) = 0, consider the coupled processes r]t and Ct with 
initial joint measure 7°°. Let 7t denote the average of 7^ over s € [0,t], and 
7 the weak limit of 74^ along some subsequence tn- Then, 7 is stationary with 
marginals /i"', and D is invariant in time under 7. Since /i"' are assumed 
to be extremal with densities ai < 02 at 0, one has ^{D) = 0. Because 
7~(D) < 7(L»), one has 7°°(-D) = 0, as desired. □ 

In the following lemma and Theorem 8, we will assume that 

there exists a one-parameter family /z", a G (0, 1), of extremal 
^^^-^ stationary -y-profile measures with /i°{?7 : 77(0) = 1} = a. 

It follows immediately from part (a) of Theorem 6 that 

jj,°',a G (0, 1), are stochastically ordered, with /i"^ < for ai < 02- 
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For ao G (0,1) and z ^ Z'^, one also has 

(49) fi°'{r]:r]{z) = l}^fi°"'{T]:r]{z) = l} as a ^ oq- 

This follows by noting that, as either a qq or a \ ao, /i" /x', where 
fj,' is stationary with fi'{r] : r/(0) = 1} = uq. When a y {a \ uq), one has 
/i' < (fi' > fj,""^); elementary reasoning similar to the path argument near 
the end of the proof of Theorem 7 will then imply /i' = /i"", and hence (49). 
(Otherwise, the coupling measure for fx' and /i"" would have discrepancies 
at some sites, but not at 0.) From Theorem 6 and (49), it follows that ^u", 
a G (0,1), are the unique extremal stationary w-profile measures. One can 
also check that, for each z G Z*^, 

fi"{ri ■.ri{z) = 1} ^ asa^O, 

(50) 

— > 1 as a ^ 1, 

since this holds at z = 0. 

Like Lemma 5, the following lemma compares stationary measures. We 
consider here a stationary coupling of /i"" with J^^q^^ i^i°^a{da), where 
and satisfy (48), and the configurations almost everywhere have dis- 
crepancies of at most one type. Then, will "dominate" and "be dom- 
inated by" Ji^Q ^-^ fi°'a{da) on invariant subsets corresponding to a G (0,ao] 
and a G (oq, 1), respectively. We employ the notation 

D~ = {{r],C) -viz) < C{z) for some z G Z'^}, 

= {{i],C)-'n{z) > Ciz) for some z G Z'^}. 

Lemma 6. Let the coupling = {r]t,Ct) be stationary, with measure 7 
having marginals /i"", ao G (0, 1), and /^q -|^^ fi°^a{da), where fj,°'° and are 
v-profile measures satisfying (48), and a is a probability measure. Also, let 
7 be concentrated on configurations with discrepancies of at most one type. 
Then, there exist invariant subsets G~ = {0, 1}^ x H~ and = {0, 1}^ x 
that partition the space, with /^q fi°^a{da) and f^^^ n°^a{da) concen- 
trated on H~ and , respectively. Moreover, 

-f{D- n G-) = -f{D+ n G+) = 0. 

Proof. The existence of invariant sets H~ and that partition {0, 1}^ 
and are concentrated on J^^ ii°'a{da) and //^^ fi°'a{da) follows, with a 
little work, from the unique representation of the stationary measures in 
terms of their extremal elements. It is easy to see that G± = {0, 1}^ X 
partition {0, 1}^'^ x {0, 1}^'' and are invariant on 7. 

We will show that ^{D~ n G~) = 0; the argument for "y{D^ n G^") is anal- 
ogous. Assume, on the contrary, that 7(1?^ H G^) > 0, and denote by 7", 
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the stationary measure obtained by conditioning on D~ n G~ . By the ex- 
tremality of the first marginal measure of j~ is /i~ = Also, by 
the extremality of /i", a G (0,ao]i the second marginal measure A~ can be 
written as 

A"= / /iV~(da) 

"'(0,ao] 

for some probability measure a~ . This implies 

fi-{rj:rj{z) = l}>X-{C:C{z) = l} 

for all z. However, since 7" is supported on , this inequality cannot hold 
for all z, which produces a contradiction. Hence, 7(-D~ H G~) = 0, as desired. 

□ 

In the proof of Theorem 8, we will employ an extension = {f]t,Ct) of 
the coupling = {r]t,Ct) introduced already, where the coordinates of fft = 
('Vf )«6(o,i) have initial measures fi" which satisfy (48), with r]Q^{z) < tj^'^iz) 
for ai < 02 and all z. It then follows that rj"^ (z) < rj^'^ (z) for all t. We denote 
by A the measure of (q, by /2 the measure of r/o and by 7 the measure of 
'^0 = (?7o)Co)- As before, we assume ffo and (q are independent. 

Theorem 8 is the main result of the section. It states that when there 
exist f -profile stationary measures /i" satisfying (42) and (48), and a v- 
profile measure A satisfying (42), then the measure of the process starting 
at A will converge to an average of the former as t goes to infinity. In one 
dimension, this statement is an immediate consequence of the convergence 
theorem for positive recurrent Markov chains. In higher dimensions, quite a 
bit more work is required to prove it. 

Theorem 8. Assume that there exist v-profile measures fx", a £ (0,1), 
that are stationary for the exclusion process and satisfy (42) and (48). Also, 
assume that (42) is satisfied by a v-profile measure A. Then, 

(51) At ^ Aoo as t ^ 00 
for some stationary v-profile measure Xoo, where 

(52) Aoo= / fi'^aida) 

and a is a probability measure. 

Proof. The argument consists of two main parts. We will first show 
that any limit of At^,, along a subsequence tk, must be of the form (52). We 
will then show that a in (52) does not depend on the subsequence, and so 
this limit must in fact hold along the entire sequence. The first part of the 
argument uses (44) of Lemma 3, (45) of Lemma 4 and Lemma 5; the second 
part uses (43) of Lemma 3, both parts of Lemma 4 and Lemma 6. 
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Characterization of limits. We will first show that any limit can be ap- 
proximated locally by a mixture of /i", a = j'/n, j = 1, . . . , n — 1, if each /i° is 
restricted to an appropriate subset of configuration space. We will then show 
that conditioning on this restriction does not change the distribution of /x" 
much. Such a limit will therefore be close to an average of the unrestricted 
Taking limits as n ^ oo will imply (52). 

We apply here the coupling that was introduced between ryf , a G (0, 1), 
and Ct, with initial measures /i" and A. By (45) of Lemma 4, for each N 
and choice of a = j /n, j = 1, . . . ,n — 1, {vjfXt) typically does not have dis- 
crepancies of both types (0,1) and (1,0) on Bn, for large enough t. Also, 
on account of (44) of Lemma 3, its analog for £'y^^{z), and (50), 

P^ii- Vf{z) > Ct{z) for ah z G L^(0)} ^ as a ^ 0, 
P'fif. r]f{z) < Ct{z) for ah z G L^(0)} ^0 as a ^ 1 

and N ^ CO, uniformly in t; since L^(0) C Bj\[, we may replace L^(0) with 
Bn in the display. It follows that, for given e > 0, large n,N and small 

as>0, 



u ''^r) 



(53) U ^^F/'") >\-e 

for large enough t, where A^^n = {j - j G [oe'i- + 1) (1 — CLs)n — 1]} and where 

Npj,n ^ ^^.^ ^,,n(^) < ^^(^^ < ^,+l,n(^) all . G S^} H (^[j ^F/' . 

That is, off of a set of probability e, the space can be partitioned into sets 
^F/'"', j G ^e.n, such that, for each j, the configuration is sandwiched 
between r]{'"' and ffl'^^'"' on Bj\i. (Here and later on, when a = j/n, we often 
write r/^'" for r/".) 

We do not know much about the behavior of rfj.'^ on ^F^'^. We will show, 
however, that rfj.;"' is almost independent of ^F^"^ for t' = t + s, where s is 
large but small relative to A^, if we consider only z G Bm, for fixed M. First 
note that, when s G [0, sq] for given SQ,\i N is chosen large enough, it follows 
from (53) and the shift invariance of the transition probabilities that 

P^i U {i^''Fr--rf^''{^)<Ct'{z)<r,l^^'\z)io,^\\zeBM}] 

(54) 

> 1 - 2e, 

for given M. That is, on ^F^'^', the configuration typically remains sand- 
wiched between ry^;" and r]^i^^'^ on Bm- Let ^Xf denote the measure ob- 
tained from Q' by restricting configurations from Z'^ to Bm, let ^^'^^j,'" 
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denote the measure obtained from rjp"' by also conditioning on ^F^"^ , and 
set ^Cj'" = Also, let •) denote the total variational distance 

between measures of configurations on Bm- (Since the space is finite, this is 
just the sum of the absolute value of the differences of the probabilities over 
all configurations.) Employing this notation, we can restate (54) as 

(55a) Y: ""^r V'-r. < *''"^A. < V?'", 

for some ^^'^Af with 

(55b) dM(*%/''^AiO<2e. 



On the other hand, for each j, one can consider the coupled process 
{v'siCs)-! where the initial measures are given by and /x-''" and the 

coordinates are initially independent with joint measure 7'. As before, it 
follows from (45) of Lemma 4 that, for given M, {r]s,Cs) typically does not 
have discrepancies of both types on Bm for large s. Moreover, since both 
measures are extremal stationary and the density of /i-''" at is greater than 
that of /x-'"^'", it follows from this and Lemma 5 that, for large s, 

P^'iiv', C) : V'siz) < C{z) for aU z G Bm} >l-e' 

for given e' > 0. Set = rjl'^ and condition on the set ^Fl'^ . For such s, it 
follows that, for e' small relative to ^cj-'"', rj^i^ typically dominates, at sites 
in Bm, a process with measure fi^~^''^ (since conditioning on ^Fj^'^ does 
not change the measure of the process starting at t/q). The same reasoning 
also shows that, after conditioning on ^F^'^, "Ht^^'^ ^^'^ large s, typically 
dominated, on Bm, by a process with measure /x-'"''^'", when e' is small 
relative to -^c^'". This domination translates into the inequalities 

(56a) V"''" < '^'''^/i^;" and < V+''" 

for some and ^^'^fl'j.J^^'^ with 

(56b) ^ci>dM(*'''^M^;"/'''^4'P < e^^ci■'"^^A/('''^/i^■,t/'^*'''^/^t^'•") < ^' ■ 

Setting e' = e/n, (55) and (56) together imply that, for given s > and 
large n, N, 

(57a) Y ^^t^'" V"''" < '-''"^V < E ^c^'"V+''" 

for some ^'^'^Xf with 

(57b) dM(^V,^'^V)<4e, 
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li s = t' — t is large relative to M, but small relative to t' . Suppose now that 
At' Aoo as t'f, — > CO, for some measure Xoo- Then, the analogue of (57a) 
holds with n = rifc ^ oo and = — > oo, but with (57b) replaced by 

dM {^^\ , *^'^At/^ ) ^ as A: ^ oo. 

Since (^^/i-'""^'", ^^^-''''^'") — > uniformly in j as n — > oo, it follows, with a 
little work, that 

*'Aoo = / V^^(rf«), 

J{0,1) 

where a is the weak limit of J2jeAe n ^'^t^^j/n along some subsequence f'^. ; 

because the measures EjgA,,„ ^4'"'5o-i)/n and EieA,,„ ^Ct'"'5(i+2)/n are 
concentrated on [a^, 1 — Og], the sequence is tight on (0, 1). {5a denotes the 
point mass at a.) Letting M— s- oo, one obtains 

Aoo = / n°'cr{da), 
J{o,i) 

which is (52). It follows from this representation for Aqo [or, alternatively, 
from (44)] that Aoo is a -y-profile measure. 

Uniqueness of a. We now show that weak limits Aoo and A'^^^ of At along 
different subsequences must, in fact, have the same measures a and a' 
in (52). Let g{ao,a), with aQ,a £ (0,1), denote the expected number of 
discrepancies of type (0, 1) on L{0) for the pair (rj, Q with marginal mea- 
sures and /i", where the joint measure is concentrated on configurations 
satisfying 

r]iz) > C{z) for all z, if oq > a, 

(58) 

< C{z) for all z, if ao < a. 
The first part of the argument will be to show that 

(59) / g{ao,a)a{da) = / g{ao,a)a' (da) 

"'(0,1) "'(0,1) 

for each ao- 

We consider the coupled process with marginal initial measures 

ao £ (Oj 1)) and A. In addition to assuming Af^ ^ Aoo as /c — > oo, with 
Aoo satisfying (52), we may also assume that the joint measures jt^ Too, 
for some 700- Let 7^°° denote the process restarted from 70° = 700, and let 
7t be the Cesaro average of 7^ over s E [0,i]. Any weak limit 7 of 7t, as 
t — > 00, is stationary, with the same marginals and /^q fi"a{da). Also, 
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by (46) of Lemma 4, 700 and, hence, 7, is concentrated on configurations 
with discrepancies of at most one type. 

Let gt(ao) denote the expected number of discrepancies of type (0, 1) on 
L(0) at time t, {chq) the corresponding expected number on L^(0), and 
g(ao) the expected number on L(0) for the measure 7. By (43) of Lemma 3, 

gt{oiQ) is nonincreasing in i, so ^00(00) *== limt^oo 5't(ao) exists. Also, by (46) 
of Lemma 4, 



Qt (o^o) — * as — > 00 
uniformly in t. Therefore, since 74^. ^ 700, it follows that 

5(00) = 5oo(ao)- 

We claim that 



(60) / g{aQ,a)a{da) = g{aQ). 

i(o,i) 

Since 5oo(«o) does not depend on the choice of t^, it will follow from the 
previous two equations that 

g{aQ,a)a{da) = / g{aQ,a)G' {da) = g^{aQ), 
(0,1) i(o,i) 

and so (59) in fact holds. 
To show (60), note that 

-g{ao) = E^[#z G L(0) : v{z) < ({z); G"] + E^[#z G L{0) : rj{z) < ({z); G+], 

where we choose G~ and as in Lemma 6. By Lemma 6, the first term 
on the right-hand side is 0. The second term on the right-hand side equals 



E 

zeL{o) 



^"{^ : ri{z) = l}a{da) — fi'^°{r] : ri{z) = l}a{a : a > oq} 

(ao,l) 



+ E^#zeL{0):v{z)>az);G+]. 
By Lemma 6, the above expectation is 0. It follows that 

5(ao) = E / fi°{rj-viz) = l}cr{da)-n'^°{i]:r]{z) = l}a{a:a> ao} 

On the other hand, it is easy to see that, because of the coupling in (58), 
I{ao 1) 9 {(^0,0:)'^ (da) is equal to the right-hand side of the last equation and 

g{aQ,a)a{da) = 0. 

(0,Qo] 

Putting these terms together, we see that (60) holds, which implies (59). 
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We will show that (59) is not possible unless a = a' . We first note that, 
for h G (0, ao), 

g{aQ — h,a) — g{a(),a) = g{aQ — h, oq) for a > oq, 

(61) 

= for a < ao — h. 

The first equality holds since uncoupled particles from the first term on 
the left-hand side arise from either uncoupled particles from the second 
term on the left-hand side or from the term on the right-hand side; the 
second inequality holds trivially since both terms on the left-hand side are 0. 
Consequently, for each ao £ (0, 1) and h G (0,ao), 

{g{ao -h,a)- g{ao,a))a{da) = g{ao - h, ao)(7{[ao, 1)), 
with the analogous equality also holding for a' . Also, note that for a G 

< g{ao -h,a)- g{aQ,a) < g{ao -h,a)< g{ao - h, ao). 
The last two displays, together with the second part of (61), imply that 

{g{ao -h,a) - g{ao, a))a{da) - / (^(ao -h,a) - g{ao,a))a' (da) 

(0,1) 7(0,1) 

< g{ao - h, ao)(cr([ao - h, 1)) - cr'([ao, 1))). 

By (59), the left-hand side of this inequality is 0, and so the right-hand 
side is nonnegative. Note that g{ao — h, ao) is at least the difference of the 
probabilities of there being a particle at for and for ^u"""'*, which by 
assumption is h. So, 

a{[ao-h,l))>a'{[ao,l)). 

Letting /i \ implies 

a{[ao, 1)) > a'{[ao, 1)) for each ao G (0, 1). 
Since the reverse inequality also holds, one has a = a' , as desired. □ 

When (6) holds, the product measures Ua^, c > 0, in Corollary 1 of The- 
orem 6, satisfy the conditions in (48); it is easy to see that they also sat- 
isfy (42). One thus has the following immediate consequence of Theorem 8. 

Corollary 3. Assume that the kernel p{-) satisfies (6), and that a 
given v-profile measure A satisfies (42). Then, 

At ^ Aoo as t ^ oo 
for some stationary v-profile measure Aqo satisfying (52). 
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8. Hydrodynamics and local limits. Hydrodynamic scaling has been ap- 
plied to a wide variety of stochastic processes, including the exclusion pro- 
cesses on Z'^, to obtain detailed information about their asymptotic behav- 
ior. In the setting of the exclusion processes, one obtains deterministic limits 
that satisfy the well-known Burgers' equation. These results extend to lo- 
cal limits for the unsealed exclusion processes, that are product measures 
with densities given by the solution of Burgers' equation there. After stating 
these results, we provide several concrete examples of such limits for exclu- 
sion processes in d = 2. We then apply hydrodynamic scaling to obtain a 
quick (though based on some nontrivial results) alternative proof of The- 
orem 4 (with the weak inequality). As in previous sections, the exclusion 
processes are assumed to be on Z'^, with transition probabilities that are 
translation invariant and irreducible; here, they are also assumed to have 
finite range. For the reader's convenience, certain basic results pertaining to 
entropy solutions of Burgers' equation are given in the Appendix. 

The time evolution of the asymmetric exclusion processes on Z with near- 
est transition probabilities p{-) has been extensively studied for product ini- 
tial measures with constant densities r and i over the positive and negative 
half-lines. (See [14] for references.) Setting p(l) = 1 — p(— 1) =p and letting 
Xt denote the corresponding measure at time t, one has, for p> 1/2, 

Ui /2 , for r < i and ^ > ^ , 

, , . Ur, for r > i and £ + r > 1, 

}^^'=\ue, for£<|and£ + r<l, 

^i^e + ^Ur, for < £ < r and ^ + r = 1. 

On the level of heuristics, the first three lines of (62) can be motivated by 
using approximations leading to the one-dimensional Burgers' equation 

(63a) ^+rn-^[u{l-u)]=0, 

with 

u{0, y) = r for y > 0, 

(63b) 

= i for y < 0, 

where m is the mean of p{-)- Under this initial data, the entropy solution 
of (63) evolves in two different ways depending on whether or not £ <r. If 
i <r, then u{t,-) is given by the shock wave that is a translate of u{0,-) 
in (63b), with 

u{t, y) = r for y > m{l — r — i)t, 

(64) 

= £ for y < m(l — r — i)t. 
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li i> r, then u{t,y) is a rarefaction wave that is continuous for t>0, with 
u{t, y) =r for y > m(l — 2r)t, 

(65) 

= £ for y < ?Ti(l - 2£)t, 

and is hnear over [m(l — 2£), m(l — 2r)]. Where M(t, •) is locally nearly con- 
stant, it is reasonable to expect A* to be close to a product measure. Sub- 
stitution of y = into (64) and (65) then gives the densities in (62). 

Hydrodynamic scaling provides a rigorous connection between entropy 
solutions of Burgers' equation and the asymptotic behavior of the exclusion 
processes in d dimensions. For exclusion processes, the relevant formulation 
of Burgers' equation in d dimensions is 

(66a) ?^ + E-^^Kl--)]=0' 

i=l * 

with 

(66b) u(0, x) = tio(x), 

with measurable uq{x) € [0,1], xG W^, and m = {mi, . . . ,mj) =J2z^Pi^)- 
Its connection with the exclusion processes is given by Theorem 9, which 
is a paraphrase of Theorems 1.3 and 7.1 in [16]. The first result (70) is the 
fundamental hydrodynamic limit; the second result (71) is a modified local 
limit. The original theorems apply to zero range processes as well. 

Theorem 9 employs the following notation. We let A", n = 1, 2, ... , denote 
product measures on Z'^, with 

(67) A"{7?:7?(z) = l}=n,,„ for^GZ^ 

where, for all r > 0, 



(68) / \u]^nx],n ~ ''J'o{x)\ dx ^ as n ^ oo, 

J|x|<r ' 

and uq{x) G [0,1] is measurable. For a cylinder function / on {0,1}^'' [i.e., 
depending only on r]{z) with \z\ <r, for some r], set 

(69) /(a) =E^" [/(??)], aG[0,l]. 
Also, set Tzf{rj) = /{tzT]), where Tzr]{z') = t]{z + z'). 

Theorem 9. Suppose that the product measures A", n = 1, 2, . . . , satisfy 
(67) and (68) for some uo{x). Then, for any finite open ball B C , i > 
and e > 0, 



(70) P 



A" 



X! Vnt{z) - / u{t,x)dx 
zenB •'^ 
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as oo, where u{t,x) is the entropy solution of (66). Moreover, for any 
cylinder function f, 



(71) 



A" 



zGnB 



nt 



B 



f{u{t, x)) dx 



as oo. 



The limit (71) states that, as n — > oo, A^^ converges locahy to a product 
measure when viewed on the hydrodynamic scale. With certain monotonicity 
restrictions on uo{x), it follows from (71) that, away from discontinuities 
of u{t,x), AJ^j in fact converges locally to a product measure without any 
averaging. This is the content of Proposition 3. 

Let Cy^e be the cone with vertex at the origin, pointing in the direction 
V ^0, and including all points in Z"^ within angle 9 >0 of v. We will require 
that the product measures A" satisfy 

(72) A"{77 : r]{z) = 1} < A"{r/ : r/(z + z') = 1} 

for all z G Z'^ and z' G C^^e, for a given choice of v and 0. That is, trans- 
lation by increases A" stochastically. After coupling the corresponding 
processes, application of Theorem 9 implies the following result with a little 
work. 



Proposition 3. Suppose that the product measures A", n = 1,2,..., 
satisfy (67), (68) and (72). Then, 

(73) r[r,x]Kt ^ ^u{t,x) as n ^ oo 

for any t> and continuity point x of u{t,-), where u{t,x) is the entropy 
solution of (66). 

We omit a proof of Proposition 3 since it is close to results contained 
in [10]. Theorem 3 there assumes Uo{x) is continuous rather than requir- 
ing (72); a corollary assumes that A" = A is fixed, with density ai in one 
octant in R'^ and another density 02 elsewhere. Both results hold for a more 
general class of particle systems that include the zero range processes as 
well. 

In order to obtain explicit limits in (73), one needs to be able to solve 
Burgers' equation (66) with its assigned initial data. Fortunately, the equa- 
tion is degenerate in the sense that its solutions are given by a {d — 1)- 
dimensional family of solutions of the one-dimensional Burgers' equation, 
along lines pointed in the direction fn = (mi, . . . , m^); this simplifies the 
computations. To state the result. Proposition 4, we choose an orthonormal 
basis vi, . . . ,Vd with vi = rh/\rh\ , let (y, W2, ■ ■ ■ , Wd) denote the coordinates 
of X with respect to vi, . . . ,Vd, and set w = {w2, ■ ■ ■ ,Wd)- (If m = 0, any 
orthonormal basis can be chosen.) We defer its proof to the Appendix. 
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Proposition 4. Suppose uq{x), x G W^, is measurable with uq{x) € 
[0,1]. Let u{t,y;w) denote the family of entropy solutions of (63a), with 
m = \rh\ and uo{y;w) = uo{x), and assume that 

(74) u{t,x)'^= u{t,y;w) 

is jointly measurable int andx. Then, u(t,x) is the entropy solution of (66a) 
with initial data uq{x). 

We provide here several applications of Propositions 3 and 4, where we 
explicitly solve for the solutions u{t,x) of (66). In each case, we examine At 
as t ^ oo, for initial product measures Ua with densities a{z) satisfying 

a{z) = r on A, 

(75) 

for given r S [0, 1] and A C Z^. Analogous examples hold in Z'^, d> 2. The 
following examples include two-dimensional analogues of the limits in (62), 
when viewed away from discontinuities of u{t, •). 

Example 1. Let A denote the half-space z^^^ > cz^"^^ for some c > 0, 
where z = {z^^\ z^"^^), and assume rh = (1,0). In (67), we can set Uz^n = r on 
A and u^^n = i on A'^, so that A" = A = i^on for all n, where a is given in (75). 
Since nA = A, we choose 

uo{x) = r on A, 

(76) 

^ ^ =i onA", 

in (68). The cone condition (72) is clearly satisfied. Consequently, by Propo- 
sition 3, 

(77) T[nx]>^nt ^ i^u{t,x) as n ^ oo, 

with u{t,x) satisfying (66). Application of Proposition 4 reduces the com- 
putation of u{t,x) to solving (63), with y = x*-^-* — cx^^^ and m = 1. The 
one-dimensional problem breaks into two cases, depending on whether or 
not £ <r. li £ <r, then n(t, y) is given by (64); otherwise, it is given by (65). 
Note that the resulting solution u{t, x) is measurable, as required for Propo- 
sition 4. 

Example 2. Let A denote the wedge A = AiCi A2, where Ai and A2 
are, respectively, the half-spaces where z^^^ > cz^"^^ and z^^^ > —cz^'^^ for 
some c > 0. As before, we assume that m = (1,0). It again follows from 
Proposition 3 that (77) holds for initial data given by (76). Application of 
Proposition 4 reduces the problem to solving (63), with y = x^^^ — c\x^'^'^\ 
and m = 1. As in Example 1, the problem breaks into two cases depending 
on whether or not i <r. Again, u{t,y) is given by either (64) or (65). 
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Example 3. Let A denote the wedge in Example 2, but assume that 
m = (0, 1) here. The hmit (77) again holds for the initial data in (76). Ap- 
plication of Proposition 4 reduces the problem to solving (63a) for m = 1 
and 

uo{y) = r on A{x^^''), 
= 1 onyl(x«)", 

where 

^(a;«) = {y:yG [-xW/c,a;(^Vc]}, 

with y = x^'^\ 

For a;(^) < 0, one trivially obtains u{t, y) = i. For x^^^ > 0, we first consider 
i > r, with the case i < r being analogous. We use standard arguments 
similar to those on pages 291-303 in [17]. One can check that the behavior 
of u{t, y) depends on whether (a) t < Ix^^'^ /c{i — r ) or (b) t > 2x^^^ / c{i — r) . 
Under (a), u{t,y) is continuous and piecewise linear except at y = {1 — i — 
r)t + x^^^ /c, with 

u{t, y)=r for y G [(1 - 2r)t - x^^^ /c, {l-£-r)t + x^^^ /c] , 

(78) =i for y G (-C30, (1 - 2£)t - x(^Vc) 

U((l-£-r)t + x(^Vc,oo). 

That is, there is a shock wave emanating from x^^^/c and a rarefaction wave 
emanating from — x^^^^ /c at i = 0, that first meet at to = 2x(^) /c{£ — r). 
Under (b), u{t,y) is continuous and piecewise linear except at 

(79) b{t) = (1 - 2e)t + 2^J2{e - r)xWt/c - x^^) /c, 
with 

u{t,y) = e - ^2{e - r)xW /ct for y = b{t), 

(80) =£ for y G (-CX), (1 - 2£)t - x^^Vc) 

U (6(t),oo). 

[We are assuming here that u{t,-) is left continuous at b{t), i.e., b{t) lies 
on the rarefaction wave.] The derivation of (79) and (80) requires some 
computation. If b{t), t>to, is the position of the discontinuity at which the 
shock and rarefaction waves meet, then 

(81) b'{t) = l-e-u{t,b{t)). 

This motion is analogous to that of the shock in (64). Moreover, because 
{t,b{t)) lies on the rarefaction wave, 

(82) ^(t,6(t)) = l(t_5(t)-xW/c), 
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which is analogous to the stretching of the wave in (65). Substitution of 

(82) into (81) gives a first-order hnear equation having the solution (79); 
together with (82), this also yields (80). The analogous formulas hold for 
the case £ < r, with —y being substituted for y in (78) and (80). 

In Theorem 4, we showed that if f G Z"^ \ {0} and p{-) has a finite mean m, 
then any v-profile measure will satisfy {rh,v) > 0. Using (70) of Theorem 9, 
we give a short alternative proof here. As elsewhere in this section, we assume 
that p{-) has finite range. (The condition v £ Z'^ from Section 5 is not needed 
here, however.) 

In order to employ (70), it is useful to be able to omit the assumption 
that the measures A" in (67) are product measures. Let the measurable map 
U :i2+ X — > P([0,1]) (the probability measures on [0,1] equipped with 
the weak topology) be a (measure-valued) weak limit, under hydrodynamic 
scaling, of a sequence of exclusion processes with initial measures A". We 
will assume any such limit satisfies the regularity condition 

(83) / dx [ \y-uo{x)\{U{t,x)){dy)^0 as t \ 

J\x\<r J[0,1] 

(a.e. in t), for all r > and a given measurable ^0(2;) S [0,1]. Kipnis and 
Landim [8] comment, on page 199, that one can substitute (83) in Theorem 9 
for the assumption that A" are product measures. (This is only stated for 
zero range processes restricted to the scaled unit torus, but it is also true in 
the present setting [11].) We will use this variant of Theorem 9 in the proof 
of Theorem 10. 



Theorem 10. Suppose that fi is a v -profile measure that is stationary 
for the exclusion process. Then, {rh,v) > 0. 

Proof. Set Uz^n = fJ-iv '■ "ni^) = 1} z G Z*^ and n G Z+. Because fj, 

is a t)-profile measure, Uz,n satisfies (68) with 

uo{x) = 1 for {x,v) > 0, 

(84) 

= for {x, v) < 0. 

Likewise, since fit = fj,, for all t, (83) trivially holds for limits U{t,x) of the 
exclusion process under hydrodynamic scaling. So, (70) is satisfied for the 
entropy solution u{t,x) of (66). Again, since fit = one must have, for all t, 

u{t,x) = 1 for {x,v) > 0, 

(85) 

= for {x, v) < 0. 

This implies {fn,v) > 0. Otherwise, as in Proposition 4, let u{t,y]w), w £ 
i?^^^, denote the entropy solutions of (63a), with m = \m\. If {m,v) < 0, 
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then (63b) is satisfied with i = 1 and r = 0, for each w. One obtains the 
solution (65), which contradicts (85). □ 

We note that one can also demonstrate Theorem 10 by using the con- 
stant sequence of product initial measures with Uz,n = 1(z,d)>0) rather than 
appealing to (83). The argument involves truncations, and so is not as quick 
as before. 

9. An example with local behavior depending on p(-). Results such as 
Theorem 9 and Proposition 3, that use hydrodynamic limits, provide sub- 
stantial insight into the asymptotic behavior of exclusion processes on Z"^. 
A more refined analysis is needed on the original scale, however, to obtain 
information about the exclusion process at the shocks of solutions of the 
corresponding Burgers' equation. For instance, the limit of At in (62), for 
the case where < i <r and i + r = 1, is given by a mixture of the product 
measures vi and Uj., since the shock can randomly be on either side of 0. 

We give here another type of example, that illustrates how the asymptotic 
behavior can depend on p(-) itself, and not just on its mean m, if d> 2. For 
this, we examine exclusion processes on with initial measure 

= 1} = 1 for zeA, 

= for z G A", 

> and z^'^^ > 0}, and with random walk 



p{z) =pi 

(87) 

= qi = 2 

for pi > 1/4, and 

p(z)=P2 for Z = 61,62, 

(88) 

= q2 = l-2p2 for z = 0. 

Under (87), rh = {2pi — |,2pi — i), and under (88), m = {p2,P2), which 
are, of course, equal for p2 = 2pi — i. Since particles can never move under 
(86) and (88), Xt = X for all t, in this case. As Theorem 11 shows, the 
behavior under (86) and (87) is quite different. 

Theorem 11. Assume that an exclusion process has transition probabil- 
ities p{-) satisfying (87) and initial measure satisfying (86). Then, for each 
j, 

(89) P^{r]:r]t{z) = 1 for some z^^^ + z^^^ <j} ^0 ast^oo. 



X{T]:r]{z) = 

(86) 

where A = {z = {z^'^\ z^^^) : z^^^ 
kernels 



for Z = 61,62, 
-Pl for 2 = -61,-62, 
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A little thought shows what, in principle, must be occurring for the process 
in Theorem 11. Particles randomly move into vertical and horizontal strips 
just to the left and below the z^^^ and z^^^ axes, along which they can drift 
up and to the right. The local departure of particles continues over time, 
emptying the plane locally and producing (89). 

The demonstration of Theorem 11 takes up the remainder of the section. 
Much of the work is contained in Lemma 7, which employs the following 
notation. We define A-' as in (86), but where A is replaced by Dj = {z : z^^^ + 
2;(2) > j}^ and set 

A, = {z : z(^) > -J, > - j, + < j} 

and 

Vj = D^j n {z : z^^'> < -j or z^^) < -j}. 

In words, Aj is the triangle with vertices at (— j, — j), (— j, 2j) and {2j, — j), 
and Vj is the half-plane below z^^^ + z^'^^ =j, excluding Aj. Inequalities 
(90) and (91) in Lemma 7 give upper bounds on the probabilities of there 
being particles in Vj and Aj for the exclusion process starting at A and A^-' . 
From now on, we consider pi > 1/4 in (87) to be fixed. 

Lemma 7. For appropriate c> 0, 

(90) P^{r]:r]t{z) = 1 for some z G Vj} < c{qi/piy^^ 
and 

(91) P^^' {r]:r]t{z) = 1 for some z G Aj} < c(gi/pi)-''/^ 
for all t. 

Proof. In order to demonstrate (90), it suffices to show the inequality 
with Vj replaced by Vj^i and by Vj^2i where 

Vj^i = Vj n {z : < ^{2)}^ y^. 2 = Vj n {z : z« > z^^)}. 

By symmetry, it suffices to do this for just Vj^i. 

Our main estimate will be a strengthened version of (44) of Lemma 3, 
with an explicit rate of decay. For this, we employ the notation 

L(z) ={z'gZ2 :(/)(') =z(2)}, 
t-{z)=E^'[#z':{z')^^'^ <z^^\ [z')^^) = z^'^'^ and77(z') = l] 

and 

l+{z)=E^'[#z':{z')^^^>z^^\ (z')(2)^^(2) andr?(z')=0]. 
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where A' is defined as in (86), but with A replaced by A' = {z:z^^^ > 0}. 
The terms L, £y and £y are special cases of those in Section 7, with d = 2 
and V = (1, 0) here. We let = f « denote the product measure with 

(92) niz) = ipi/qif' 

[where tt{z) = a{z)/{l — a{z))]. 

Using Theorem 2, it is easy to check that z/ is stationary. Clearly, and 
A' are f-homogeneous measures for v = (1,0). They also satisfy both parts 
of (42). Setting ii = v and A = A', the assumptions of Lemma 3 are therefore 
satisfied. Consequently, by (44), 

l-,[z)^{) asz^^^^-oo 

uniformly in t. Using the explicit form of 7r(-) in (92), this can be strength- 
ened by setting e = ci{qi/pi)"', with appropriate ci > 0, in the first two 
displays of the demonstration of (44). Combining these two inequalities, as 
before, implies that 

(93) i-,{z)<2ci{qi/pir 

for z^^^ = —2n and all t. Since Ac A' , it follows that 

(94) £l{z)<2ci{qi/pir 

for z^^^ = —2n and all t. 

The region Vj^i is the union of the horizontal line segments (— oo, r(2;(^))], 
e z, with 

r{z^^^) = -3-l for G [-j, j], 
Together with (94), this implies 

oo 

i?^[#zGUj,i:r/(z) = l]<(2j + l)2ci((7iM)t-''/'l+2 2ci(gi/pi)["/2] 

n=j+l 

< c(gi/pi)^/3 

for large enough c. This demonstrates (90). 
The reasoning for (91) is similar. Here, we set 

(95) L{z) = {z' G Z2 : (z')^!) - {z')^^) = - z(2)}. 

The expectations (z) and (z) are defined as before, but with the inner 
equality in (95) replacing {z')^"^^ = z^'^\ The vector v = (1,1) replaces (1,0), 
and we let i^^ = Va denote the product measure with 

(96) 7r(z) = (piM)^*''+^''^-^-. 
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The measures are stationary; and A-' are u-homogeneous and satisfy 
both parts of (42). Setting = v^^ and A = A^-', it follows as before that, for 
appropriate ci > 0, 

£;,,(z)<2ci(giMr 

for z^^^ + z^^) = 2j — 2n and all t. The region Aj is the union of 2j + 1 line 
segments pointing in the direction v with right endpoint z^^^ + z^^^ = j in 
each case. Setting n = [j/2] in each case implies that 

E>^'' [#z e Aj : riiz) = 1] < (2j + l)2ci(giM)b'/2] < c(gi/pi)jV3 
for large enough c. This demonstrates (91). □ 

The proof of Theorem 11 applies Lemma 7, together with a coupling 
argument that employs the translation invariance of D2j under (—1,1). The 
basic idea is that since 

{z:z(i) + z(2)<j-} = y^.UAj, 

(89) will follow from (90) and (91), if the initial measure X^^ in (91) can 
be replaced by A for large t. The difference A — D2j of the sets of sites 
initially occupied by the corresponding processes is finite. So, it is reasonable 
to expect that the particles originally there will eventually dissipate away 
from the finite region Aj as t — > oo, perhaps by moving more or less in the 
direction (—1,1). 

The coupling we will use, like that after (50), is an extension of the 
coupling = ivtXt) to a coupling = (vtXt) with multiple coordinates. 
Here, ff = {rj^, . . . ,77^""^), with = [19j"*/e], for a given e > 0. The coor- 
dinate processes are assumed to have deterministic initial measures /i^-''", 
n = 0, 1, . . . , — 1, and A^-' . The measures are defined as before; ^^^'^ 
are defined as in (86), but with A replaced by D2j^n = U A2j^^, where 

A^. „ = {z : z(^) > 0, > 0, z^^^ + < j} + (-2jn, 2jn) 

is the translation, by (— 2jn, 2jn), of the triangle with vertices at (0, 0), (0, j) 
and (j, 0). Let 7 denote the joint initial measure of the coupling. 

We note that /x^-^'" is the translate of fi'^^'^ by ^; = (— 4jn,4jn), and so 
^t'^'^ is also the translate of /if"'''' by z; since A C D2jfl, ^'''^ stochastically 
dominates A. Also, since D2j C ^2j,n, for all j and n, the discrepancies for 
a pair {rj'^Xt) are always of type (1,0). Their number remains constant over 
time, which is 2j^ + j, the number of sites in A^- ^. The distribution of such 
discrepancies, at a time t, is the same as for n = 0, up to a translation by z. 
We also set 



Aj-„ = Aj + (-4jn,4jn), 
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and note that A"^ n A"^ = for ni i^ni. 

Proof of Theorem 11. Set 

Ft = if- V?{z) = 1, Ct{z) = for some z G Aj>}, 
for n = 0, 1, . . . , — 1. We will show that, for fixed e > 0, 

N-l 

(97) J2P'i^n<eN 

n=0 

for large enough t. It follows from the above comments on discrepancies that 
the probabilities of the sets in (97) are all the same, and so P^{F^) < e. 
Consequently, 

{-q : rjt{z) = 1 for some z G Aj} — P {rj : rit{z) = 1 for some z G Aj} < e, 

for such t. Choosing e < c{qi/piy^^ , it follows from this, (91) and A C i^2j,0) 
that, for given j and large enough t, 

P^{r]:r]t{z) = 1 for some z G Aj} < 2c{qi/piy^^. 

Along with (90), this implies that 

P^{ri : r]t{z) = 1 for some z^^) + z^^^ < j} < 3c(gi/pi)^'/^, 

which implies (89) and, hence, Theorem 11. 
In order to demonstrate (97), we set 

MO = #{n:ieFn. 

We claim that, on Xt > 9j^, there must exist ni,n2 and 21,2:2) with Zi G 
Aj^m , so that 

(98) r^^^(z,)=v?'iz2) = l and r?^ (z^) = (^1) = 0- 

To show this, let n'l, . . . ,n'^p denote the first indices for which ^ G 
Then, 



3j 



2 



U{^:V(^) = 1, 0(^)=0} 



i=l 

So, there is an index ni, with 



<3j2|A,| <9/. 



(99) rfl'{zi) = l and r/"'(2i) = 0(21) = 

for some 21 G Aj^n^ and all i. For some 722 G {n'^, . . . ^n'^p}, 

{2 G Aj- „2 : '?rn-z) = 1 and Ct{z) = 0} = 0, 
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again since |Aj| < 3f. One can therefore choose Z2 G ^j,n2 so that 

r?r(^2) = l and ^J^^ (^2) = Ci(^2) = 0. 

Together with (99), this imphes (98). 

When the event in (98) occurs, the pair (?/"^,f/"^) has discrepancies of 
opposite types at zi and Z2- Since Zi S Aj^„., one has \z2 — zi\ < 4\/2jN. In 
the same manner as in the proof of Lemma 4, one can employ Lemma 3.1 
of [1] to obtain a uniform lower bound 5 > on the probability of these or 
other discrepancies of opposite type meeting over {t,t + 1]. Also, there are 
not more than Sj'^N sites at which ^0 has discrepancies. Therefore, the event 
> 9j^ cannot occur with probability at least e/2 at more than 6fN/6£ 
times spaced at least distance 1 apart. It follows that, for given e > and 
sufficiently large T^, 

(100) P^Xt > 9f} < e/2 

fovt>Te. 
One has 

J2 P\Ft) = E'^lXt] = E^[Xt-Xt > 9f] + E'^iXuXt < 9j% 

n=0 

Since Xt < N, this is, by (100) and our choice of N, 

< leN + 9f < eN, 
for small e. This implies (97), and completes the proof of Theorem 11. □ 

10. Open problems. An important open problem in the context of this 
paper is to determine completely the set of extremal stationary measures 
for an irreducible translation invariant exclusion process on Z'^,d>l, with 
J2x \ x\p{x) < CO. A complete characterization is out of reach at the present 
time. It is therefore useful to isolate particular parts of this problem that 
are deserving of attention. Here are some (related) statements that we think 
are probably true, but cannot prove. [Some may require higher moment 
assumptions on p{-)-] 

1. All extremal inhomogeneous stationary measures are t;-homogeneous for 
some V ^ satisfying {J2x^P(.^)^'^) > 0. 

2. All extremal inhomogeneous stationary measures are 't;-profile for some 
V & R'^ satisfying {J^x^Pi^)^'^) > 0. (We are less sure about this state- 
ment.) 

3. For each v ^ satisfying {J2xXp{x),v) > 0, there is a continuous one 
parameter family of inhomogeneous extremal stationary u-profile mea- 
sures. 
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4- If J2x xp{x) = 0, then the extremal stationary measures are exactly {z^p, < 
p < 1}. This is an old problem — it appears as open problem ^5 on 
page 416 in [13]. 

5. For every choice of p{-) with nonzero mean vector, there exists an extremal 
stationary measure that is not a product measure. [We do not know that 
this is the case for any p{-) if d > 2.] 

In one dimension, ^1 above is trivial, ^2 is true, ^3 is true with extra as- 
sumptions, except that the parameter is discrete rather than continuous, ^^4 
is true, and #5 is true under extra assumptions. (For the extra assumptions 
required, see the discussion of the one-dimensional system in Section 1.) 

Problems #1, and ^,^3 are interesting and open if Z'^ is replaced 
by a "ladder" of the form Z x {1,...,A^}. A possible analogue of the v- 
homogeneous and f-profile properties in this context is the following: One 
can view {1, . . . , N} as a cycle, and say that a measure is rotationally homo- 
geneous if it is invariant under "rotations" of Z x {1, . . . , N}, and is profile if 
it is rotationally homogeneous and has density tending to 1 in one direction 
and to in the other. Problem ^4 can be handled in the case of a ladder 
by a small modification of the one-dimensional proof. It is likely that ^5 for 
the ladder can be done using the approach in [2]. 

It may be useful to recall the known characterization of the class of all 
extremal stationary measures for continuous time independent particle sys- 
tems on Z*^ with transition kernel p(-), since this might shed some light on 
what to expect for the exclusion process. By Theorem 4.12 in [12], the ex- 
tremal stationary measures are exactly the Poisson fields on {0, 1,2, . . 
in which the number of particles at x has mean m{x), where 771(3;) > for all 
X and J2x m{x)p{y — x) = m{y) for all y. By the Choquet-Deny theorem, the 
extremal functions m satisfying these properties are the pure exponentials 

m(x) = e<^'^'\ 

for V such that J2x e~^^''"^pix) = 1. Nonconstant m(-) of this type can only ex- 
ist if J2x xp{x) 7^ 0, and then the corresponding v must satisfy {J2x xp{x),v) > 
0. Of course, nonextremal m's can correspond to extremal stationary Poisson 
fields. 

Note that the Poisson fields corresponding to the pure exponentials are 
u-homogeneous (and even satisfy a property analogous to being f -profile), 
while those corresponding to sums of exponentials are not. In one dimen- 
sion, there are no extremal stationary measures for the exclusion process 
corresponding to sums of exponentials — they are all either homogeneous or 
v-profile. Perhaps this is true in higher dimensions as well. This has a bearing 
on problem ^1 above. 

In the context of the convergence results of Section 7, it is not clear to 
us what assumptions on p{-) and v are, in general, necessary in order for 
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the conditions (42) and (48) of Theorem 8 to hold. It seems reasonable to 
guess that the existence of a one-parameter family /x" of extremal stationary 
v-profile measures, as in (48), is in some sense generic. 

The tail condition (42) should presumably hold for the ?;-profile measures 
of a large class of p(-). A related question arose in [1] in d = 1. As pointed out 
there, presently there are no known examples of extremal stationary profile 
measures that are not blocking measures. Our condition (42) corresponds 
to (1.3) in that paper, which immediately implies that a measure is a block- 
ing measure. It is also not known what conditions on p(-) are necessary for 
all stationary blocking measures to satisfy (1.3); it is suggested there that a 
third moment assumption on p{-) might be the correct condition. 

APPENDIX 

In Section 8, we examined the entropy solutions of the d-dimensional 
Burgers' equation in (66). The equation is a specific case of the scalar con- 
servation law 

(Ala) |^ + ^„.^;,(„)=„, 

1 = 1 ' 

with fi G C^{R), and having bounded initial data 
(Alb) u{0,x)=uo{x), xeR^. 

We briefiy review here basic existence, uniqueness and stability results for 
the entropy solutions of (Al). We then give the quick proof of Proposition 4. 
More detail can be found on conservation laws in one dimension in Chap- 
ters 15 and 16 in [17] and Chapter 3 in [4], and on conservation laws in 
d-dimensions in [9] and in Appendix 2 in [8]. 

Classical solutions to (Al) need not exist for all t, even when the initial 
data uq[x) is smooth. One therefore typically works with weak solutions 
of (Al). That is, for every function g : R^^ ^ ^ R oi class C]^^{R+ x i?'^), 
u : i?+ X i?'^ — > is required to satisfy 

u{t,x)—g{t,x) + ^mifi{u{t,x))—g{t,x) dxdt 



(A2) 



uo{x)g{0, x) dx = 0, 



where fi and uq are as in (Al). (Cj^^ denotes the continuous functions with 
compact support and one continuous derivative in both time and space.) 

Solutions u{t,x) of (A2) can be modified on sets of measure on x 
R'^ without affecting the left-hand side of (A2), and so are certainly not 
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pointwise unique. Because of the possibility of shocks, they also need not 
be a.e. unique. (See, e.g.. Chapter 15 in [17] for examples.) The physically 
meaningful solution of (A2) is given by its (a.e.) unique entropy solution. 
The following entropy condition is due to Kruzkov [9], and states that, for 
every positive function g of class (ii+ x R'^) and every a > 0, u{t,x) 
satisfies 



, \\u{t,x) - a\-^g(t,x) 
JB'^ ot 



(A3a) 



d (9 1 

+ yi"^i|/i(^(*;3;)) - fi{a)\—g{t,x) \dxdt> 0, 



and for all r > 0, 



(A3b) / |u(t,x) — uo(x)| ^ ast\0 

(a.e. in t), where ft and no are chosen as in (Al). A solution u{t,x) of (A2) 
satisfying (A3) is said to be an entropy solution. (Kruzkov [9] includes a 
more general family of equations in place of (A2).) 

Kruzkov [9] proved the existence and uniqueness of entropy solutions: 



Theorem Al. For every bounded measurable uo{x), there exists a unique 
entropy solution of (A2). 

The above uniqueness is a consequence of an stability result in [9] that 
compares entropy solutions with different initial data. 

The theory for solutions of (A2) in one dimension preceded that for gen- 
eral d (see, e.g., [15]). When, in addition, /" < 0, there is a version of the 
entropy solution satisfying 

(A4) u{t,y)-u{t,y + a)<Ca/t, 

for a,t> 0, where C > depends only on / and ||no(y)||oo (see, e.g., page 
266 in [17]). In particular, only increasing jumps can occur as one moves 
from left to right, and for a given t, u{t, •) is locally of bounded variation. In 
this setting, one can employ (A4) as the entropy condition in place of (A3). 

We conclude with the proof of Proposition 4 of Section 8. The argument 
is a straightforward application of Fubini's theorem. 



Proof of Proposition 4. We need to show that (A2) and (A3) hold 
for fi{u) = f{u) = u{l — u) and all g{t,x), if u{t,x) satisfies (74). Using the 
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coordinates of the orthonormal basis introduced before Proposition 4 and 
Fubini's theorem, we can rewrite the left side of (A2) as 



where g{t,y;w) = g{t,x). Since u{t,y;w) is an entropy solution of (63a) and 
g{-,-;w) G C]^^{R+ x R) for each w, the quantity inside the braces always 
equals and, hence so does the entire integral. So, u{t,x) satisfies (A2). The 
argument for (A3) is the same. □ 

We note that the joint measurability of u{t,x) in t and x was needed in 
order for (A2) to make sense and to apply Fubini's theorem. By altering 
each one-dimensional solution u{t,y;w) at a single point {tw,yw), one can 
produce nonmeasurable u{t,x). So, one cannot arbitrarily choose the version 
of u{t, y; w) for each w. We also point out that the above proof does not hold 
when fi depends on i, since one would need to factor the functions outside 
of the summation in order to obtain the analogue of (A5). 
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